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Section  1 


INTRODUCTION 


Ionospheric  disturbances  cause  severe  disruptions  of  radar  and  communi- 
cations systems,  and  high  background  rates  in  optical  defense  systems.  The 
most  troublesome  disturbances  affecting  high  frequency  wave  propagation  are 
the  scintillations,  which  can  occur  at  any  time,  and  are  expected  to  be 
especially  severe  after  a high-altitude  nuclear  explosion.  An  extensive 
program  is  now  underway  to  find  out  how  ionospheric  disturbances  affect 
radiowave  propagation  through  the  ionosphere.  The  program  comprises  experi- 
mental studies,  such  as  WIDEBAND,  on  naturally  occurring  disturbances  in  the 
ionosphere,  and  theoretical  studies  of  the  scintillations  and  of  the  ionospheric 
conditions  which  lead  to  this  occurrence.  This  report  addresses  several 
aspects  of  the  theoretical  investigation  of  ionospheric  irregularities.  It 
is  mainly  concerned  with  the  effects  of  particles  and  fields  in  the  upper 
ionosphere  — a problem  which  must  be  solved  if  the  underlying  causes  of 
irregularities  are  to  be  understood.  The  emphasis  in  the  first  sections 
is  on  processes  in  the  natural  ionosphere.  The  last  section  treats  an  effect 
that  may  be  expected  to  occur  after  a high-altitude  nuclear  detonation. 

A new  method  was  recently  described  for  calculating  the  distribution  of 
electrons  leaking  out  of  the  permanently  trapped  radiation  belts  (Ref.  1,  2). 

The  results  are  important  for  their  implications  to  magneto spheric  physics 
and  the  coupling  of  the  ionosphere  to  the  magnetosphere.  An  important 
by-product  of  the  calculation  is  the  energy  deposited  in  the  upper  ionosphere 
by  electrons  precipitating  in  mio 'latitude  regions  (L  < 6).  The  distribution 


of  enemy  dep_;.i.tion  can  be  calculated  quite  accurately  because  the-  method 
incorporates  the  AURORA  code,  which  has  been  carefully  developed  and  tested 
■ver  a period  of  many  years  (Ref.  3,  k,  1).  It  was  therefore  suggested  that 
results  from  this  program  could  be  of  great  value  to  tie  interpretation  of 
data  from  WIDEBAND  and  other  experiments  directed  at  the  assessment  of 
radio  propagation  through  the  ionosphere.  The  fundamental  advantage  is  that 
the  precipitation  of  electrons  from  the  trapped  particle  belts,  and  the  re- 
sultant energy  deposition,  are  more  understandable  and  readily  predictable 
than  in  the  auroral  regions.  Eventually  it  might  be  possible  to  predict 
the  particle  precipitation  contribution  to  ionospheric  heating,  even  without 
we  11 -coordinated  satellite  observations  --  something  that  will  probably  never 
be  accomplished  for  the  auroral  regions.  The  advantages  have  been  somewhat 
offset  by  the  lack  of  mid-latitude  observational  data.  The  sparcity  of 
particle  precipitation  data  has  precluded  the  calibration  of  diffusion  and 
precipitation  rates.  The  mid-latitude  WIDEBAND  data  were  limited  mainly  to 
one  station  that  only  operated  for  a portion  of  a four  month  interval, 
though  we  intend  to  look  into  the  possibility  of  using  low  elevation  data 
from  the  WIDEBAND  receivers  at  high  latitudes  to  investigate  the  northern 
edge  of  the  mid-latitude  region. 

Sections  2 and  3 describe  a continuation  of  previous  work  (Ref.  1) 
and  first  attempts  at  interpreting  WIDEBAND  data  using  our  results.  Section  2.1 
is  a discussion  of  the  basic  theory  behind  the  calculations.  It  has  heretofore 
been  assumed  that  the  conventional  bounce -averaged  diffusion  equation  might  be 
adequate  to  treat  the  electron  precipitation  problem.  We  shall  show  how  our 
approach  is  related  to  the  bounce -average  method,  and  why  a bounce -averaged 
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treatment  can  never  yield  satisfactory  results.  Section  2.2  describes  further 
efforts  to  calibrate  the  electron  diffusion  rates  by  observations  of  trapped 
and  precipitating  particles.  The  results  have  noL  been  very  encouraging, 
mainly  because  of  instrumental  difficulties.  However,  some  of  t:he  difficulties 
might  be  overcome  by  a new  method,  outlined  in  Scotian  2.3,  that  takes  advan- 
tage of  the  variation  in  electron  fluxes  due  to  their  longitudinal  drift 
motion.  The  following  section,  3,  is  concerned  with  a preliminary  analysis 
of  WIDEBAND  data.  1'nfortunatel.y,  those  data,  and  the  supportive  data  on 
magnetic  variations,  ionospheric  densities,  and  trapped  particles  did  not  be- 
come available  soon  enough  to  allow  a detailed  statistical  analysis  for  possible 
correlations.  There  were,  however,  a sufficient  number  of  interesting  scintilla- 
tion events  at  mid-latitudes  to  lend  support  to  our  suggestion  that  studies 
of  ionospheric  activity  and  particle  precipitation  could  profitably  be  carried 
out  in  mid- latitudes. 

Recent  satellite  measurements  of  the  angular  and  energy  distributions 
of  ions  and  electrons  by  Shelley  et  al.  (Ref.  5)  and  Mcllwain  (Ref.  6)  and 
elec;  ‘Lc  fields  of  Mozer  et  al.  (Ref.  7)  provide  improved  information  on 
phenomena  affecting  the  coupling  of  the  ionosphere  and  magnetosphere  and  on 
the  acceleration  of  ion  from  the  upper  ionosphere  to  energies  in  the  keV  range. 

The  reason  these  results  are  of  interest  here  is  that  they  imply  the  widespread 
and  frequent  occurrence  of  plasma  Instabilities  that  lead  to  anomalous  resistivity 
and  electric  "double  layers,"  which,  in  turn,  cause  large  electric  fields 
aligned  along  the  magnetic  field.  These  phenomena  may  lead  us  to  modify  our  en- 
tire understanding  of  the  exchange  of  particles  between  the  ionosphere  and 
magnetosphere, and  of  the  motion  of  nuclear  debris  after  a high  altitude  nuclear 
explosion. 
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Section  4 describes  an  interpretation  of  some  ion  and  electron  dis- 
tributions observed  >y  Shelley  et  al.  (Ref.  5).  The  events  selected  were 
those  wherein  high  fluxes  of  keV-type  h"*"  and  0"*"  ions  were  observed  moving 
away  from  the  atmosphere  and  highly  directed  along  the  magnetic  field.  In 
the  analysis,  an  electric  field  was  assumed  extending  below  and  above  the 
satellite  along  the  magnetic  field.  In  the  calculations  described  here  the 
distribution  of  ions  and  electrons  was  followed  downward  from  the  observation 
point,  through  the  electric  field  region,  into  the  atmosphere,  and  back  up 
to  the  observation  point.  The  AURORA  code  was  used  to  compute  the  flux  of 
backscattered  electrons.  The  computational  model  was  much  simpler  than 
that  described  in  Section  2 — ignoring  the  pitch-angle  diffusion  due  to 
wave-particle  interactions  at  high  altitudes.  The  observations  could  be 
explained  quite  convincingly  by  passage  through  a potential  drop  of  as  much 
as  4.5  kV.  The  data  also  Indicated  that  high  fluxes  of  keV-type  electrons 
were  trapped  between  the  converging  magnetic  field  below  and  the  electric 
field  above. 

One  of  the  prime  candidates  for  a mechanism  that  produces  ionospheric 
irregularities  is  the  Far ley-Buneman  instability  (Ref.  8)  that  is  driven 
by  the  kinetic  energy  in  the  streaming  motion  of  Lhe  plasma.  If  the  stream- 
ing velocity  is  greater  than  the  thermal  velocities  the  instability  grews 
rapidly.  Obviously,  a high-altitude  nuclear  explosion,  with  the  attendant 
plasma  motions,  is  a likely  candidate  for  plasma  instabilities.  The  atmospheric 
heave  drives  plasma  across  the  magnetic  field,  inducing  a large  electric 
field.  This  electric  field  polarizes  the  ionospheric,  and  the  polarization 
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charges  propagate  to  high  altitudes  along  the  magnetic  field.  The  resulting 


electrostatic  field  extends  to  large  areas  of  the  ionosphere,  far  beyond  the 
region  principally  disturbed  by  the  burst.  The  Far ley-Buneman  instability 
can  be  excited  if  this  electric  field  in  the  effective  E-region  of  the  iono- 
sphere induces  a strong  drift  motion  of  the  electrons  with  respect  to 
the  ions  --  at  a rate  higher  than  the  ion  thermal  velocity. 


Section  5 describes  a test  of  the  Far ley-Buneman  instability  criterion, 
using  some  MRHYDE  and  MICE  computer  data  on  upper  atmosphere  conditions  follow- 
ing two  H.A.N.E.  events.  Unfortunately,  all  of  the  data  required  for  a complete 
analysis  were  not  available.  Nevertheless,  the  analysis  indicated  that  the 
electrostatic  field  arising  from  the  atmospheric  heave  is  indeed  very  large 
and  that  magnetic  field-aligned  plasma  Irregularities  due  to  this  instability 
may  be  expected  to  extend  to  several  thousand  kilometers  from  the  burst. 
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Section  2 


DIFFUSION  OF  TRAPPED  PARTICLES 
AND  THE  DISTRIBUTION  OF  PRECIPITATING  PARTICLES 

2.1  THE  PITCH-ANCLE  DIFFUSION  EQUATION 

The  motions  of  trapped  particles  iri  the  magnetosphere  can  be  studied  by 
a statistical  treatment  that  results  in  a Fokker- Planck  diffusion  equation. 

By  considering  times  longer  than  the  particles'  gyroperiods  one  of  the  six 
degrees  of  freedom  can  be  eliminated.  Another  variable  can  be  eliminated  by 
either  assuming  uniformity  of  the  distribution  in  longitude,  or  by  acknowledging 
that  the  longitudinal  drift  motion  is  predictable,  and  linear  with  time.  That 
leaves  four  parameters  (plus  time)  to  describe  the  particles'  distribution  — 
one  advantageous  choice  is  L-shell,  latitude  or  position  on  a field  line, 
energy,  E,  and  pitch  angle,  a.  We  may  also  choose  to  ignore  the  L-parameter 
and  concentrate  on  the  physical  processes  that  involve  only  the  remaining  three 
parameters.  We  do  not  mean  to  imply  that  L-shell  — or  radial  — diffusion  has 
no  significant  effect  on  the  distribution  of  precipitating  particles.  Radial 
diffusion  is  important  and  should  eventually  be  included  in  the  treatment  of 
particle  losses  (Ref.  9).  The  remainder  of  this  section  is  devoted  to  practical 
methods  for  solving  the  Fokker -Planck  equation,  primarily  by  judicious 
elimination  of  one  or  more  variables. 
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2.1.1  The  Bounce -Averaged  Diffusion  Equation 

Stably  trapped  particles  interacting  with  l->w- level  waves,  or 
with  i tier  particles  i>b--y  a di  tYusii.n-type  oquati  <k:  (Ret'.  10,11). 

dr 

*-  V cos  a constant  ad.iriv. 


3t 


(1) 


^ l.D...  a in  a ^-t]  - [energy  loss  term]  + [sources]  - r losses] 


— : ■< — ’ A J .aj.ii  Lx  ^ — 

sin  a act  ~ act  oa 

wlior"  f(t ,s ,K,a)  is  a phase  spac:  number  — density  distribution  function,  and 

D , (b,s,E,oj  is  a pitch  auric  diffusion  coefficient.  The  remaining  independent 
oi(y 

variables  are  time,  t,  and  distance  along  a field  line,  s,  (usually  measured 
from  the  equator);  V is  the  velocity.  The  energy  loss  term  is  important  for 
collisions  with  particles  in  the  atmosphere  (Ref. 10),  but  can  generally  be 
ignored  for  interactions  of  electrons  with  waves  at  high  attitudes 
Kef.  11,  12).  Most  of  the  important  physical  processes  are  described  by 

the  first  term  on  the  right  side,  so  we  will  concentrate  on  that  term,  and 
ignore  the  remaining  three  terms  (usually  small)  except  when  the  source  and 
loss  terms  arc  needed  for  logical  consistency. 

Equation  (l)  is  a local  diffusion  equation.  The  second  term,  however, 
must  be  differentiated  with  the  first  adiabatic  invariant  hold  constant.  A 
suitable  replacement  for  the  local  pitch  angle  is  any  variable  that  depends 

only  on  the  adiabatic  invariant.  A suitable  adiabatic  invariant  variable  is, 

\ 

if  energy  remains  constant; 


S0  2 
u = 1_  bR  sin  “ 


(2) 


where  (s  = 0 ) is  the  minimum  magnetic  field  strength  - or  equatorial  field 
and  B(a)  is  the  local  field.  For  free,  non-interacting  charged  particles, 
u would  be  i tirrns  the  cosine  squared  of  the  pitch  angle,  a^,  at  the  equator; 
in  the  present  case  there  is  no  simple  functional  relation  between  u and  cr^. 
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Thu  diffusi  equation  can  be  written 
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Th’1  conventional  method,  for  treating  this  equation  had  been  to  average  ever  the 
lesgth  of  tin;  trajectories  between  the;  mirror  points,  s and  a,..  The  element 
■ t"  distance  along  a particle  trajectory  is  ds/cos  a.  A simple  integration 
gives 


b - m 


1 ^ 
T 5u 


ds 

cos  a 


1 i 9 
J T 5u 


where  the  subscripts  1 and  2 indicate  values  at  the  conjugate  mirror  points; 
< > denotes  an  average  over  ds/cos  »,  and  T is  the  "quarter  bounce  period" 
integral 


ds 

cos  a 


(5) 


The  constant  R has  the?  dimensions  of  a Length;  in  a dipole  field  it  is  just 
the  radial  distance  to  the  equatorial  crossing.  In  a (nearly)  symmetric  dipole 
field  a simple  empirical  approximation  holds  for  T (Ref.  13). 


T(u)  **  1 . i8017 j 


•6jyf>93  (l-u)3/8  ' 


(6) 
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Ir  RqutiVifii  (•'*)  th--  commutativity  ft’  the  operations  ^ arid  < > is 
hi 


and  r—  remain  bounded  at  the  mirror  points.  The  separation  of 
era-  do 


insured  i d ' P( 

S ,, 

Oi 

I and  in  the  average  on  the  right  side  would  further  simplify  the 
HU  Ou 

■quatio!!,  though  this  roqu  1 res  an  assumption  that  is  not  always  warranted 
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uu 

*8] 

(7) 

[Note  that  an  average  over  ds/cos  a is  the  only  way  to  remove  the  hS/hs  term, 
though  the  original  .justification  for  the  averaging  procedure  did  not  take 
this  into  account  (R'f.10)].  In  the  loss  cone,  a < a , the  assumption  of 
uniformity  of  f along  the  field  line  cannot  be  valid.  In  the  limiting  case 
of  a totally  absorbing  atmosphere,  where  no  particles  are  scattered  backward 
into  the  upward  direction,  the  df/ds  term  becomes  Vf^/T,  The  resulting  equation, 
which  has  been  frequently  used  for  the  less  cone,  is 


|i-  ss  | I <D  > T 

ot  T ou  I uu 


f <D  > T jpl  - — 
L uu  Ou  J t 


where  is  the  bounce  period  (4RqT/V).  Somewhat  different  derivation  have 

been  used  elsewhere  to  justify  Equation  (8).  An  alternate  method  is  to 

assume  that  the  particles  arc  lost  within  one-fourth  bounce  period,  on  the 

average.  The  loss  turn  is  then  4f/T  (Ref.  14,  15,  16,  17).. 

B 


It  is  relatively  easy  to  demonstrate  that  the  bounce  averaging  procedure 
runs  into  difficulties  when  the  variability  of  f with  s is  tak-  n into  account. 
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The  bounce  avc-r-  cd  diffusion  equation  is  now 


9 < f > V , v 
“~3t  + T f i 

= ^ I-  [<  D > T « » 

T du  L uu  du 


[<  D > T « »1 

f uu  ou  J 

i r<B  >T(L«i_2>  ,1  «:».t  <f>\] 

T ou  L uu  ou  L /J 

where  the  symbols  « » denote  a new  Kind  of  average 

s , B s , B, 

« X » - f “ X ~ cos  a ds  / T ~ cos  a as 

vO  •it) 

S,  3, 


(9) 

(10) 


(il) 


(1-) 


Thus  an  attempt  to  remove  one  integration  variable,  s,  has  onLy  resulted  in  the 
appearance  of  a new  dependent  variable,  « f 5^,  ou  the  right  side  of  the 
diffusion  equation.  It  may  be  concluded  that,  except  for  special  cases,  the 
bounce-averaged  diffusion  equation  is  not  valid  for  the  treatment  of  distribu- 
tions in  and.  near  the  loss  cone.  Bounce-averaging  is  only  justified  for  slow 
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!'  '' permanently"  trapped  p-artic  1 es  with  large  pitch  angles. 


u i 1’  f lit:  i < ’ 1 1 


A rough  criterion  car:  be  formulated  for  the  range  nf  validity  of 
Equation  (7).  If  there  is  a cutoff  equatorial  pitch  angle,  o\  . below  which  all 
the  ti’aipid  particles  interact  strongly  with  the  atmosphere  within  half  a 
bounce  period,  the  distribution  can  be  symmetric  in  s only  for  pitch  angles 
a greater  than  n + Ao  , where 

a ■“  0 ) 


Aa 


J T . < D_ 


Of  y 

O 'vJ 
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(ii) 
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v B 
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"aa  F 


/ t_  < D . ^ ^.s-?  > 


cos  a 


Values  of  < D > are  available  from  several  sources  (e.g.  Ref.  19,20);  normally 
a a 

rj  o 

D -k  2 -l 

< at  a > is  expected  to  be  less  than  10  rad  sec  . Table  I shows  the 

o o 

numerical  factors  that  give  Ay  from  < D >.  Except  for  the  case  of  strong 

o O'  or  r 

o o 

diffusion,  when  < D > becomes  comparable  with  l/Tg,  Equation  (7)  is  satis- 

CJ  o 

factory  for  most  of  the  trapped  particles.  But  near  and  within  the  loss  cone. 

0fQ  £ or  , the  pitch  angle  diffusion  equation  must  be  solved  explicitly  for 
(at  least)  three  independent  variables. 


To  see  what  can  be  accomplished  with  Equations  (/) and (8) , consider  the 
idealized  case  where 


D . = D (sin  ot  )** 

sin  a , sin  a o o' 

o o 


(14) 
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TABLE  I 


Range  o£  Influence  of  Atmosphere  on 
Permanently  Trapped  Particle  Distribution 


E ev 

Att  (electrons) 

o 

(protons) 

io3 

l.bSyJh  < D > 

11.1^/l  < D ^ 

10* 

.953  " 

6.26  " 

105 

.568 

3.52  " 

I06 

.434  " 

1.98  " 

io7 

.420  " 

1.12  " 

th.-it 


and  til' 


maintains  the  trap-p”'!  particles  is 


S (sin  a ) 1 


(1^) 


wh-T1.'  D and  B are  ••constants . The  diffusion  equation  for  a stati  nary  distri- 


bute (crf/an.,  0)  is,  after  rep. I ncing  sin  or  by 


1 i S 

trapped:  '■>=■=•  - 


tb  * ^ 


loss  cone:  0 


f ",  _ q*-  .L  dfl  „ p-ll  . 

LDoV  3xrso*  ix>\ 

1 1 3 fn  - <i+1  3 r 1 2f  | 

Tn  X *X  lDo  B*  oxj  t 1 x _ Xc 


(16a) 

(lbb) 


To  simplify  the  solutions  let  be  a constant,  equal  to  the  bounce  period 


at  > . The  solutions  of  Equation  (16a) are  straightforward,  and  give  for 


trapped  particles 


f = - r 

o 
S 


_ i Vp+^‘q  + c V_q  + C 1 d / o 

D (P+2)(p^-q)  X °1  X C2  ! ^ * 


f = - 


o 


~~2  XP+"  *-  CL  In  X + C£l  q = 0 


(17a) 

(17b) 


o (p +2)‘ 

where  the  C's  are  constants.  The  loss  cone  equation,  (lob) , is  Bessel's  equation 


and  gives  solutions  in  terms  of  Bessel  functions,  ln,of  imaginary  argument 


f*C3*“l/2 


2 / 2y2~q\1//2  I / 

q/2-q  2-q  \ i q 


(18a) 


f = C3/X 


| q = 2 


(18b) 


It  is  necessary  to  match  the  functions  and  their  first  derivatives  at  XQj 
the  remaining  constant  can  be  eliminated  by  letting  f = 1 at  x - 0.  The 
solutions  for  q ^ 2 are 


u •-.) 


1/L= 


Tiles-'  results  are  fairly  insensitive-  to  the  strength  of  the  cure-.  a great 
simplification,  therefore,  occurs  in  the  limit  as  S becomes  zero;  the  loss 
cone  solutions  ar.- 
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(20b) 


q i o • 


These  are  similar  to  the  solutions  for  the  same  case  found  by  Theodor idis  and 
and  Paolini  (Ref.  14)  except  that  they  normalized  their  solutions  somewhat 
differently. 


Several  solutions  for  the  weak  source,  q = 0 case  are  shown  in  Figure  1. 
The  figure  -was  drawn  for  50  kev  protons,  which  are  not  reflected  by  the 
atmosphere,  on  h = 4.  The  reason  for  considering  protons,  is  that  the 
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LOG  (FLUX) 


2.1.2 


Solution  of  the  Three-Dimensional  Diffusion  Equation 


Th<-  integration  of  Equation (3) poses  some  novel  difficulties,  mainly 
because  of  tin-  distortion  of  the  domain  of  integration.  The  numerical 
to'shiii'jU'-s  were  described  previously  (Ref.  l),  so  the-  following  discussion  is 
i1  ■ voted  mainly  t--  fundam-nbal  quest i >ns. 


The  time  integration  of  the  pitch  angle:  diffusion  equation  is  not 
always  interesting,  and  can  be  avoided  by  assuming  either  a steady  state,  or 
a lowest  eigenmode  decay  with 

f - g (u,s)  e"t/T  (^l) 

The  characteristic  decay  period  is  t . The  differential,  equation  for  g is 


3, 


o a K 

v ^ L ~ 
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5u 


«] 


C&) 


uu  V i-(i-u)b/b 

where  v is  ti'.e  dimensionless  distance  from  the  equatorial  plane,  S/R  . The 
domain  of  integration  (for  d dipole  field)  is  shown  in  Figure  2.  The  extreme 
limits  of  v are  ± T(0)  •-  ± 1. 300173.  The  domain  of  integration  is  bounded 
on  the  right  by  the  curve  v = v^(u),  where  the  subscript  m denotes  that  a 
parameter  is  evaluated  at  the  mirror  point.  The  figure  shows  only  half  of  the 
domain;  the  earth's  field  is  sufficiently  symmetric  that  the  Lower  half  can  be 
considered  a reflection  of  the  upper  half.  The  particles  generally  enter  the 
atmosphere  at  different  values  of  v,  say  -v  and  +v  , in  '.he  tw  ■ opposite 
hemispheres,  so  the  integration  is  not  really  symmetric.  Several  horizontal 


dashed  lines  have  been  drawn  to  indicate  for  different  L shells  the  v values 

c 
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the  diffusion 


wh‘‘iv:  n1.:no;-,phn ri scattering  and  energy  loan  must  bn  brought  into 
equation . 


Tin-  transformation  from  local  pitch  angle , a,  to  adiabatic  invariant 
variable  u,  has  altered  the  domain  <>f  integration  in  another  way  that  is  not 
immediately  apparent  in  Figure  2.  Everywhere  except  at  the  equator  (v  - 0)  the 
distribution  in  u has  been  split  into  two  separate  parts:  a group  of  particles 

going  in  one  direction  with  a % arc  sin  J S (l-u)  and  another  group  going  in 

'O  O 

the  opposite  direction  with  a J?  tt  - arc  sin  J — (l-u).  The  choice  of  u ^ cog‘‘  a 

o v B o 

o 

rather  than  ~ cos  has  caused  the  separate  regions  to  lie  on  top  of  one 
another,  though  they  are  actually  joined  only  at  the  point  v = 0,  u = 1.  A simple 
way  around  the  difficulty  is  to  label  the  two  distributions  with  superscripts 
UP  and  DOWN,  and  perform  the  v integration  in  each  direction. 


v 


If  f(o')  is  to  be  continuous  at  a - ~ 


v» 

m 


UP,  V 
g (u,vm) 


DOWN  , » 

8 («.»J 


it  follows  that,  on  the  line 


(23) 


If  the  domain  were  symmetric  about  v = 0,  only  one  g would  need  to  be  considered 
with 

g(u,Vm)  = g(u,  - vm)  (24) 


This  condition  leads  immediately  to  the  conclusion  that  g cannot  be  separated 
into  a function  of  v multiplied  by  a function  of  u (as  f was  separated  to 
give  g and  a function  of  t).  If  g were  separable  the  v eigenfunctions 
would  have  the  form  exp(-Cv);  which  is  incompatible  with  Equation  (24) except 
in  the  trivial  case  C * 0 (or  f independent  of  s). 
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It  ala'i  appears  that,  if  3,-/3 v is  to  remain  bounded  at  v - v the  doriva- 

tive  3g/3u  must  vanish  at  the  mirror  points.  For  most  of  the  interesting  pitch 

angle  diffusion  mechanisms,  the  diffusion  coefficient  approaches  zero  so  fast 

at  « = n/,_  that  the  two  distributions,  ^ and  , can  be  considered  fully 

independent.  Even  if  D uoes  not  approach  zero  at  least  as  fast  as  cos  a,  the 

ofy 

derivative  condition  follows  from  the  differential  equation,  22. 


The  boundary  condition  at  u = 1 is 


g = #>  “ 0 

w yJ 

This  does  not  imply  that  3g/3u  is  zero,  though  this  type  of  boundary  condition 
on  g(u)  is  fairly  easy  to  incorporate  in  a numerical  finite  difference  integration. 


The  most  interesting  boundary  is  at  v = v^;  here  the  integration  switches 
over  from  generally  slow  pitch  angle  diffusion  at  nearly  constant  energy  to  rapid 
diffusion  and  energy  loss  due  to  collisions  with  atmospheric  atoms  and  molecules. 
The  time  scales  in  the  two  v regions  are  so  vastly  different  that  it  is  imprac- 
tical to  solve  the  combined  diffusion  equations  in  a routine  fashion,  integrating 
over  both  regions  simultaneously  until  convergence  is  attained.  The  character- 
istic times  for  diffusion  in  the  magnetosphere  are  hours  to  days;  the  character- 
istic times  for  diffusion  in  the  atmosphere  are  less  than  the  particles'  bounce 
periods  - of  the  order  seconds  or  less.  It  is  just  this  disparity  that  renders 
liable  to  suspicion  any  attempt  to  treat  the  distinct  physical  processes  on  an 
equal  footing.  For  example  Spjeldvik  and  Thome  (Ref.  21  ) attempted  to  treat 
the  atmospheric  interactions  by  introducing  a diffusion  coefficient  in  the  bounce- 
averaged  diffusion  equation  that  was  very  large  in  the  loss  cone  region.  Apart 
from  the  difficulties  of  constructing  a numerical  solution  across  the  boundary 





where  some  parameters  change  by  orders  of  magnitude,  their  results  do  not 
appear  entirely  reasonable.  At  the  transition  between  the  two  modes  of  diffu- 
sion they  show  a sudden,  steep  drop,  followed  by  a nearly  isotropic  distribu- 
tion within  the  loss  e ns.  Yet  one  would  expect  the  downgoing  flux  at  the  top 
of  the  atmosphere  to  show  a nearly  exponential  fall-off  at  the  transition,  with 
a characteristic  scale  given  by  Equation  (13) . In  Figure  3 [similar  to  Figure  3 
of  Ref.  1 except  that  the  no-backscatter  ,-ase  shows  the  new  results  of  Equa- 
tion (19) lare  compared  to  the  Spjeldvik  and  Thorne  results  and  ou>’  new  results. 
For  a wave  field  of  30  my  the  predicted  diffusion  coefficient  for  50  kev 
electrons  is  approximately 

D > *>  10_l4  sec"1  (26) 

o u 

This  gives  a AQQ  of  ab».ut  .5°  at  the  equator.  The  distribution  at  300  kev 
altitude  should  therefore  be  spread  out  over  nearly  18°  near  the  "cutoff"  at 
75°.  The  results  of  Spjeldvik  and  Thorne  fall  off  much  faster  than  the  expected 
exp  (of/l8°)  rate;  our  results  are  consistent  with  the  expected  fall-off  rate, 
and  are  nearly  indistinguishable  from  the  bounce-averaged  results  where  those 
results  are  expected  to  be  applicable  (a  £ 6o°) , 

A few  words  of  caution  are  in  order  here.  While  scattering  within  the 
atmosphere  may  be  treated  in  the  numerical  integration  as  a boundary  condition 
at  v * v , it  is  not  generally  possible  to  construct  a simple  boundary  condi- 
tion there.  Consider  the  simple  case  where  no  particles  are  assumed  to  be 
reflected,  or 

gUP  (u,vc)  = 0 (27) 
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Figure  3.  Comparison  of  bounce -averaged  diffusion  calculations  with  new  results.  The  curve 
marked  SPJELDVIK  and  THORNE  is  from  Ref.  21.  The  curve  labeled  NO  BACKSCATTER  is 
a plot  of  Equation  (19). 


This  u-'i.ditloi.  is  iiiuomjjatibl*:  with  Equation (23) unl« :s s i-  aluo  7, urn  at 

v = v , which  in  turn  -iuni-is  tho  possibility  of  particle  loss  in  the  atmosphere. 

Tli"  difficulty  is  removed  if  the  particles  which  would  mirror  a small,  but 

finite,  distance  beyond  v are  returned  after  boiritf  slightly  attenuated. 

A simple  prescription  based  on  this  observation  can  be  used  in  most  situations 

whore  the  physical  processes  are  understood.  Often  the  "penetration  depth" 

is  small  enough  that  the  mathematical  difficulties  only  affect  one  point 

v (u. ) = v it:  an  array  of  discrete  points, 
mi'  c * 

Ideally,  the  integration  would  be  joined  across  the  boundary  at  v , to  a 
detailed  solution  of  the  Fokker  Planck  equation  describing  scattering  in  the 
atmosphere.  The  AURORA  code  generally  suffices  for  that  purpose  except  that  the 
change  in  the  pitch  angle  due  tc  magnetic  reflection  was  incorporated  in  the  code 
only  in  an  approximate  fashion.  To  the  best  of  our  knowledge,  none  of  the  atmospheric 
scattering  and  energy  deposition  codes  have  explicitly  taken  into  account  the 
precautions  discussed  above.  As  a result  it  has  taken  some  experimentation  to 
establish  the  optimum  altitude  for  the  transition  to  atmospheric  scattering. 

In  our  previous  calculations  an  altitude  of  150  to  180  kin  was  chosen.  At  higher 
altitudes  (above  200  km),  errors  quickly  build  up,  even  in  the  parts  of  the  distri- 
bution that  are  not  strongly  affected  by  the  atmosphere.  This  means  that  it  might 
not  always  be  feasible  to  obtain  reliable  results  for  electron  energies  below  1 keV. 

The  case  of  zero  reflection  is  very  much  like  what  actually  happens  to 
trapped  protons.  Protons  entering  the  atmosphere  are  rapidly  degraded  to  thermal 
energies  before  suffering  appreciable  deflection  (Ref. 10).  The  steady  state  differ- 
ential equation  (22)  has  been  solved  for  a range  of  values  of  D up  to  the  limit  of 
strong  diffusion.  The  most  readily  accessible  published  values  of  the  diffusion 
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The  simplest  assumption 


coefficient  (e.g.  Ref.  19)  are  for  the  average  <Dao«0>* 
about  the  form  of  D^,  is  that  it  is  a function  only  of  or  . It  is  curious  that 
this  assumption  leads  to  values  of  < > that  are  just  slightly  less  than 

one-half  D over  most  of  the  range  of  a . The  conversion  formula  is,  from 
Equations  6 and  9 


D ^ > - 

a a 

o o 


_ 2 
T cos  a 


aa 


(28) 


.332  + .169  tan^or  (sin'^o^  - l) 


_ 3/j4 

1 + .46 j sm  ' or. 


act 


Trie*  diffusion  coefficient  < Dq  > was  also  assumed  independent  of  or^ 

o o 


Some  sample  results  of  calculations  for  50  keV  protons  on  L = 4 are 
shown  in  Figures  4 through  7-  These  results  should  be  compared  with  Figure  1, 
which  treats  the  same  case  by  the  bounce  averaging  method.  However,  it  is 
now  possible  to  calculate  the  distribution  at  any  point  along  the  field  line. 
Even  if  the  bounce-averaged  diffusion  equation  were  fully  applicable  it  would 
only  give  a sort  of  "average"  distribution.  In  the  cases  shown  we  have  plotted 
the  distribution  at  the  equator  and  the  distribution  at  the  top  of  the  atmos- 
phere (in  terms  of  an  adiabatic  invariant  parameter  Bm/Bo  - l/l-u). 


At  small  values  of  B /Bo  there  is  a discrepancy  that  is  simply  a con- 
sequence of  stopping  the  calculation  before  the  solutions  for  particles  with 
large  pitch  angles  had  fully  converged.  In  this  region  the  bounce-averaged 
equation  gives  better  results  with  less  demands  on  computer  time.  The  new 
results  beyond  Bm/BQ  ~ 50  are  satisfactory,  and  should  be  accepted  as  the 
correct  solutions  against  which  others  should  be  judged.  It  is  peculiar  that 
for  moderate  diffusion  rates  the  averaged  diffusion  equation  results  would  agree 
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LOG  (FLUX) 
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Figure  6. 


101  201 
B (MIRROR) /B0 

Loss-cone  distributions  at  50  keV  protons  on  L=4;  < D > 

10"6sec"L  The  two  curves  shown  are  for  the  aoao 

distribution  at  the  equator,  and  for  the  distribution  going 
downward  at  the  top  of  the  atmosphere. 
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The  fact  that  the  b ounce-averaged  proton  distribution  si  are  s:.  near  tc  the 
correct  solutions  in  the  loss's  >ne  might  load  one-  to  suppose  that  the  bounce 
averaged  equation  could  work  as  well  for  electrons,  giver;  a satisfactory 
heruistic  loss  term.  There  are  several  factors  that  work  against  this  supposi- 
tion. First,  there  is  the  expectation  that  the  distribution  at  very  small 
pitch  angles  should  be  almost  entirely  due  tc  the  atmospheric  scatter i>**  processes, 
which  are  exceedingly  difficult  to  incorporate  in  a simple  model.  Second,  the 
energy  loss  of  electrons  in  the  atmosphere  must  be  considered.  It  is  unlikely 
that  calculations  at  a single  energy  such  as  reported  by  Spjeldvik  and  Thorne, 

(Ref.  21)  could  ever  yield  more  than  a crude  order-of -magnitude  estimate  of  the 
total  precipitation  rate  (though  their  distributions  are  questionable,  compari- 
sons with  uur  new  results  indicate  that  both  methods  predict  approximately  the 
same  number  of  particles  in  the  loss  cone). 


The  major  conclusions  that  can  be  drawn  at  this  point  are  that  our  basic 
method  is  sound,  and  that  the  only  significant  uncertainties  are  in  the 
predictions  of  backseattered  electron  distributions  by  the  AURORA  code.  The 
close  agreement  of  the  proton  distributions  with  results  of  bounce-averaged 
calculations  indicates  that  the  numerical  integration  method  is  giving  correct 
results.  That  the  bounce  averaged  formulation  is  good  over  such  a large  range 
of  pitch  angles  means  that  the  seemingly  crude  assumption  for  the  form  of  D 

Q(Q 
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2.2 


TESTS  OF  THE  PITCH-ANCLE  DIFFUSION  THEORY 


A careful  observation  of  the  distribution  vf  electrons  within  the  loss 
:.=ti"  luis  great  potential  as  a moans  of  do  tor  mining  both  diffusion  rates  and 
atm.-.sph'  rio  .^aeksuattor  rates.  Inspection  of  Figures  h through  7 and  ■)  of 
Iv-feroii-.";  1,  (a Is  sc-'  Ref.  confirms  that  the  loo 3 c-no  distribution  is 
iii>i"ed  v-.-v-i  sensitive  t • the  diffusion  rat':.  But  the'  important  measurable 
parameters  are  net  neesr-ari  ly  those  that  simple  intuition  might  lead  or  2 
to  believe  u ;i:ful  for  this  purpose.  For  example,  total  flux  in  the  center  of 
the  loss  m sue  (or  ratio  of  or  ~ 0°  flux  to  ‘K>°  flux  at  low  altitude)  is  nearly 
independent  of  diffusion  rate,  except  in  strong  diffusion.  The  distribution 
near  the  center  of  the  less  cone  is  m"St  valuable  for  the  information  it 
pr  .'Vides  about  bachseatter  rates.  Ideally  the  ratio  <»f  downgoing  flux  to 
uppoing  flux  at  intermediate  pitch  angles  in  the  loss  cone  would  give  directly 
the  instantaneous  value  of  the  diffusion  ooeff ieient,  a quantity  that  ha j not 
yet  been  found  by  any  other  method. 

Unfortunately  most  of  the  loss  cone  data  we  have  examined  do  not  appear 
to  be  reliable.  The  data  are  very  sparse,  and  of  uncertain  quality.  It  is 
likely  that  most  of  the  data  were  contaminated  by  cosmic  ray  induced  breinsstrahlung: 
Sp.jeldvik  (Ref. 22)  has  noted  similar  difficulties.  The  only  case  we  found  in 
which  bremsstrahlung  contamination  of  the  loss  cone  counting  rates  could  be 
ruled  out  was  the  one  reported  in  Figure  9 of  Reference  2' . The  data  for  that 
case  were  collected  during  a time  of  intense  magnetic  activity  (K^  = 8-); 
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s<-  it  iu  whether  nth  or  procosues;  could  have  affected  the  trapped 

electrons • 

We  ha vo  pro at  hopes  that  the  instrumental  contamination  problems  may 
ha vo  boon  i-»vi.-v>!-{i»o  in  tin:  most  mount  satellite  •ix|juriments;  ami  we  are* 
avia i tin i'  the  new  data.  But,  in  view  ul'  the  pravt:  difficulties , wo  wish  to 
uuppost,  in  the  foLlewiru'  soot  Lon,  another  method  that  rnlpht  be  used  to  derive 
diffusion  cneffii;  rents. 


39 


2.3  LONGITUDINAL  DRIFT  AND  THE  FILLING  OF  THE  LOGS  CONE 

An  important  effect  that  has  been  ignored  in  the  preceding  discussion 
is  the  variation  of  the  size  of  the  loss  cone  with  geomagnetic  longitude. 

The  earth's  magnetic  field  is  only  very  roughly  approximated  by  a simple 
dipole  field.  For  most  purposes  the  important  deviations  from  the  simple 
model  can  be  accommodated  in  an  off-center  dipole  field  model  (Ref.  23). 

The  cutoff,  , then  varies  with  azimuthal  position.  Figure  8 shows  the 
100  km  altitude  trace  plotted  against  geomagnetic  longitude  for  two  L shells. 
The  particle's  azimuthal  drift  allows  them  to  diffuse  beyond  or  and 

partially  fill  the  so-called  "quasi  trapping"  regions  where  the  cutoff  is 
reduced.  A diagonal  line  is  sketched  in  Figure  8 to  indicate  how  the  "edge" 
of  the  distribution  advances  as  thy  particles  drift  (in  this  sample,  electrons 
drifting  eastward) . 

If  the  diffusion  coefficient  is  large  enough,  the  diffusing  particles 

can  advance  into  the  gap  between  a and  a as  fast  as  a recedes.  The 

e c-max  c c 

maximum  change  in  a ^ with  longitude  i is  of  the  order 


|f  ~ -35/L2  . 


(29) 


The  mean  atQ  distance  a particle  can  diffuse  during  a complete  bounce  period 
ia  given  by  Equation  (13)*  The  diffusion  coefficient  near  needed  to  keep 
the  gap  filled  is  therefore 
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< It  > . Equation  (3:))  however,  refers  only  to  the  necessary  diffusion 
<y  a 


coefficient  that  will  allow  a small  number  of  particles  to  leak  out  of  the 
trappier  region  at  every  longitude . As  the  diffusion  coefficient  is  increased 
it  will  gradually  approach  a value  that  follows  f to  ad.just  continuously  to 
the  changing  cutoff,  with  a nearly  uniform  leakage  rate  independent  of  longitude. 
This  latter  limiting  value  is  much  less  than  the  strong  diffusion  limit, 

< D > ~ ^ , but  as  much  as  5 decades  higher  than  the  values  of  Table  II. 


Figures  9 and  10  show  results  of  sample  calculations  using 

Equation  (7)  for  the  diffusion  of  electrons  (with  an  assumed  < D > = 

°o  o 

constant)  with  no  allowance  for  baokseatter.  Similar  results  have  appeared 
previously  (Ref.  24).  For  low  values  of  the  diffusion  coefficient,  the  tail 
of  the  distribution  fills  in  nearly  as  cxp( -const  ~_);  the  local  cutoff,  B /B  , 

i5  CO 

O 

is  generally  out  of  sight  at  the  lower  right  of  the  figure.  When  the  diffusion 

-5-1  . 

coefficient  becomes  large  enough, in  this  case  10  sec  , the  pitch  angle  distri- 
bution is  similar  at  every  longitude,  except  that  the  local  cutoff  determines  the 
location  of  the  shoulder  of  the  distribution. 


The  filling  and  depletion  of  the  quasi-trapping  jss  cone  has  been 
observed  for  trapped  electrons  and  for  trapped  protons.  Figure  H shows  some 
results  of  Imhof  (Ref.  25)  for  a low  altitude  satellite  whose  east-west  motion 
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TABLE  II 

Minimum  Pitch  Angle  Diffusion  Coefficient  That  Allows 
Some  Filling  of  the  Loss  Cone  at  All  Longitudes 


(ev) 

D^m  (electrons) 

Dlim(pr°tons) 

io3 

4. 3 — /L 

1.9"— /L 

io4 

1.4““/L 

5 . 9*— /L 

io5 

4. 2**“/L 

1.9"-  /L 

io6 

1.3'-  /L 

5.9'~  /L 

io7 

7.4'-  /L 

1.9'-  /L 
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OMNIDIRECTIONAL  FLUX 


b/b0 


I Figure  9.  Solutions  of  the  averaged  diffusion  equation  for  various 

longitudes  (relative  to  the  SA  anomaly,  denoted  by  the 
label  on  each  curve).  This  case  is  for  a small  diffusion 
coefficient  so  the  bottom  edge  of  the  particles  remains 
nearly  fixed  at  one  B value. 

i 


i 
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Figure  10.  Solution  of  the  averaged  diffusion  equation  for  various 
longitudes  (relative  to  the  SA  anomaly,  denoted  by  the 
label  on  each  curve).  This  case  is  for  a diffusion 
coefficient  large  enough  that  the  edge  of  the  distribution 
follows  a constant  altitude  trace. 
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Figure  11.  Filling  of  the  "quasi  trapping"  electron  loss  cone  observed  by  a satellite  at 
constant  altitude.  The  solid  curves  are  results  of  a diffusion  calculation 
by  Imhof  (Ref.  25). 


was  adequate  to  exhibit  the  effect.  We  have  done  a simple  calculation  using 
Equation  8 with  the  loss  term  multiplied  by  1 minus  a "reflection  coefficient", 
R.  A plausible  fit  to  the  data  is  shown  in  Figure  12.  The  backscatter  of 
electrons  does  tend  to  mask  the  eastward  drift  effect;  but  since  the  quasi 
trapping  phenomenon  affects  mainly  particles  mirroring  above  the  atmosphere, 
there  is  hope  that  the  effect  may  be  of  value  in  determining  diffusion 
coefficients . 


How  then  might  the  atmospheric  backscatter  be  treated  simultaneously 
with  longitudinal  drift, while  avoiding  the  rigors  of  a diffusion  equation 
with  four  independent  variables,  longitude,  latitude  pitch  angle, and  energy? 
Again,  it  should  be  feasible  to  take  advantage  of  the  fact  that  the  various 
pehnomena  have  quite  different  time  scales.  If  the  cutoff  boundary,  ® , 
does  not  move  too  rapidly,  then  the  steady  state  equation,  22,  should  give 
valid  results  for  the  atmospheric  backscatter.  The  averaged  diffusion  equa- 
tion should  then  be  adequate  to  treat  the  major  part  of  the  trapped  distribu- 
tion, plus  the  quasi  trapping  region.  This,  f i^st-approximation  could  be  used 
as  a starting  point  for  several  iterations  of  Equation  22  to  obtain  detailed 
solutions,  for  all  longitudes.  We  expect  to  test  this  method  soon  to  find  out 
whether  it  is  practical,  and  indeed  results  in  great  savings  of  computer  time 


RANGE  OF  MEASURED 
FLUXES  FOR  h^;_  >20 


Section  3 


OBSERVATIONS  OF  SCINTILLATIONS  IN  MID  LATITUDES 


Ionospheric  scintillations  have  been  observed  at  all  latitudes, 
though  research  on  their  effects  has  been  concentrated  in  the  equatorial 
and  auroral  regions.  Recently,  there  has  perhaps  been  a tendency  to  emphasize 
the  equatorial  regions  because  those  regions  are  most  predictable.  However, 
it  might  be  argued  that  the  mid-  and  high- latitude  ionospheric  disturbances 
more  closely  approximate  the  effects  induced  by  high-altitude  nuclear  explo- 
sions. We  have  therefore  directed  our  efforts  at  a search  for  correlations 
in  mid  latitudes  between  scintillation  activity  and  energetic  particle 
precipitation.  The  mechanisms  underlying  such  correlations  are  poorly  under- 
stood, so  our  work  has  been  mainly  explorational. 

A receiver  station  in  the  WIDEBAND  program  was  operated  at  Stanford, 
California  from  May  25,  1976  to  mid-September  1976.  The  purpose  of  the  Stanford 
station  was  mainly  to  check  out  procedures,  and  verify  the  operation  of  the 
system.  Operation  of  the  Stanford  station  was  therefore  somewhat  irregular, 
covering  only  about  three  days  of  each  week.  Nonetheless,  an  appreciable 
amount  of  data  was  acquired  showing  high  degrees  of  scintillation  activity 
(Ref.  26).  e.  J.  Fremouw  of  SRI  graciously  allowed  us  access  to  the  data, 
and  with  the  assistance  of  S.  Matthews  (also  of  SRI),  we  were  also  to  construct 
our  own  summaries  of  the  data. 

Table  nilists  the  times  of  the  observations  and  the  estimated  levels 
of  activity.  The  scintillation  activity  levels  from  A-quiet  - to  B-intense  - 
pertain  to  the  combined  amplitude  displacement,  Sv>  and  the  phase  duration, 
ap,  of  the  VHF  signals.  Our  assessments  were  made  independently,  and  may  differ 


slightly  from  those  arrived  at  by  the  SRI  group  (Ref.  26).  The  levels  A 
to  B correspond  roughly  to  the  designations  very  quiet  and  very  active, 
respectively,  in  Ref.  26,  though  we  have  tried  to  break  down  the  individual 
pass  records  by  the  elevation  angles.  Some  low  elevation  intense  scintilla- 
tions should  be  discounted  because  they  reflect  merely  a large  total  electron 
content  over  a long  path. 

We  have  also  listed  available  data  on  K.^  indices,  solar  sector  crossings, 
and  the  character  of  the  total  electron  content  (TEC)  as  observed  by  the 
WIDEBAND  receiver. 

The  first  thing  to  be  noticed  is  that  moderate-to-intense  activity  in 
mid  latitudes  can  occur  at  all  times  of  the  day,  and  with  no  obvious  pattern. 
(The  early  morning  passes  are  labelled  with  a N because  they  started  in  the 
north,  and  the  afternoon  passes  starting  in  the  south  are  labelled  with  a £.) 
Ihere  may  be  a correlation  between  the  north-to-south  gradients  of  the  iono- 
spheric electron  content,  but  the  predictive  value  of  this  correlation  is  not 
clear.  The  absolute  levels  of  the  total  electron  content  have  not  yet  been 
established,  either  by  comparison  with  (very  limited)  ionospheric  sounder 
data,  or  by  absolute  calibration  of  the  WIDEBAND  data  (Ref.  26),  if  the 
correlation  is  established  and  confirmed,  it  is  not  clear  how  it  would  be 
exploited  to  obtain  local  predictions  of  scintillation  levels.  Perhaps  new 
ionospheric  sounding  stations  might  be  established  at  critical  sites.  In 
any  event  the  electron  density  and  temperature  are  the  most  likely  ionospheric 
parameters  that  might  reflect  the  deposition  of  energy  from  the  trapped 
particle  belts.  The  three-fold  correlation  of  scintillation  activity, 
Ionospheric  conditions,  and  precipitation  of  trapped  electrons  is  a topic 
which  calls  for  an  extended  investigation. 
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We  had  initially  hoped  to  find  correlations  between  the  global 

indicators  --  such  as  K --  and  the  mid  latitude  scintillations.  Table  III 

P 

lists  the  Kp  index  and  the  character  of  the  interplanetary  magnetic  field. 

Only  the  sum,  £K  , of  the  8 three-hour  K 's  is  shown  because  activity  was 
P P 

very  low  during  the  period.  Values  less  than  20  for  indicate  quiet 

days,  and  only  values  greater  than  30  should  be  considered  significant. 

The  interplanetary  field  during  a half-day  is  characterized  by  the  direction 
of  the  field  near  the  earth;  T stands  for  toward  and  A for  away.  An 
inspection  of  the  most  readily  available  indicators,  and  the  solar 
sector  structure  reveals  no  such  correlations.  Analysis  of  a larger  body 
of  data  might  reveal  low-level  correlations,  but  such  correlations  would 
be  of  little  value  to  an  operational  system. 

We  are  just  now  beginning  to  receive  new  data  on  electron  precipita- 
tion rates,  though  we  have  not  found  any  that  are  directly  applicable  to 
the  period  in  question. 

It  should  be  kept  In  mind  that  the  WIDEBAND  experiments  were  initiated 
at  the  bottom  of  the  solar  minimum.  Magnetospheric  and  ionospheric  activity 
have  been  very  low,  but  are  expected  to  rise  rapidly  through  the  next  several 
years.  That  the  anticipated  correlations  failed  to  materialize  is  not  surprising 
because  no  great  solar  events  occurred  during  the  period  while  the  Stanford 
receiving  station  was  in  operation.  It  has,  however,  been  demonstrated 
that  high  levels  of  scintillation  activity  do  occur  frequently  in  mid-latitude 
regions;  frequently  enough  that  the  study  of  mid-latitude  scintillations  is 
a worthy  endeavor. 
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Section  4 


COUPLING  OF  IONOSPHERE  AND  MAGNETOSPHERE 

4.1  ANOMALOUS  PHENOMENA  AFFECTING  COUPLING 

Many  processes  of  importance  to  defense  problems  are  sensitive  to 
ttie  coupling  of  the  ionosphere  and  the  magnetosphere.  This  coupling  is 
severely  altered  when  plasma  instabilities  form  anomalous  resistivity , or 
electric  double  layers, in  the  upper  ionosphere  that  impede  the  flow  of 
thermal,  plasma  currents  along  the  magnetic  field.  The  resulting  electric 
fields,  with  components  along  the  magnetic  field,  preclude  the  neutraliza- 
tion of  large-scale  electric  fields  produced  in  the  magnetosphere  by  high 
altitude  nuclear  explosions.  Such  obstructions  thereby  affect  (i)  the  growth 
rate  of  plasma  irregularities  all  along  the  magnetic  tubes  containing  the 
debris  and  "atmospheric-heave"  plasma,  (ii)  the  dynamical  motion  of  the  debris 
tube,  hence  the  eventual  distributions  of  the  debris  and  the  trapped  electrons, 
and  (iii)  the  energy  deposition  in  the  upper  atmosphere  and  the  attendant 
interference  with  radio-wavs  propagation  and  optical/IR  systems.  Furthermore, 
the  obstructions  affect  natural  processes  in  the  auroral  zones, which  are 
examined  to  understand  phenomena  produced  by  high-altitude  nuclear  explosions. 

Anomalous  resistivity  and  electric  double  layers  are  due  to  plasma 
instabilities  that  inhibit  the  mobility  of  thermal  electrons  along  the  magnetic 
field.  Although  much  theoretical  effort  has  been  devoted  to  the  investigation 
of  the  pertinent  instabilities  (see  Perkins,  Ref.  27;  Kindel  and  Kennel, 

Ref.  28;  Block,  Refs.  29,30;  Carlqvist  and  Bostroo,  Ref.  31;  Papadopoulos 
and  Coffey,  Ref.  32;  Swift,  Ref.  33;  and  Swift  et  al.,  Ref.  34),  many 
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uncertainties  remain  (see  Papadopoulos , Ref.  35,  for  a recent  review  of 
this  subject). 

Much  experimental  evidence,  however,  is  now  available  on  the  electric 
fields  resulting  from  the  instabilities.  The  fields  were  inferred  initially 
from  observations  of "monoenergetic"  electrons  in  the  auroral  2ones  (see,  e.g., 
Hoffman  and  Evans,  Ref.  36;  Evans,  et  al..  Ref.  37;  Paschman,  et  al.,  Ref. 

38;  Maehlum  and  Mocctue,  Ref.  39;  Arnoldy,  et  al.,  Ref.  40;  and  Basqued,  et. 
al..  Ref.  41).  More  recently  they  have  been  inferred  from  observations  of 
energetic  ions  generally  moving  along  the  magnetic  field  away  from  the  upper 
ionosphere  (see  Hultquist  et  al.,  Ref.  42;  Reme  and  Bosqued,  Ref.  43; 

McXlwain,  Ref.  6;  and  Shelley  et  al..  Ref.  5).  In  addition,  the  electric 
fields  have  been  inferred  from  shaped-charge  barium  releases  (Wescott  et  al.. 
Ref.  44;  Haerendel  et  al.,  Ref.  45)  and  direct  measurements  (Mozer  et  al., 

Ref.  7;  Mozer,  Ref.  46).  The  data  indicate  that  the  electric  potentials 
extend  to  6 kV  and  higher,  and  are  located  within  a few  thousand  kilometers 
from  the  earth  (Evans,  Ref.  47;  Lampton,  Ref.  48;  Winiecki,  Ref.  49;  Arnoldy, 
Ref.  50;  and  Haerendel  et  al..  Ref.  45),  Moreover,  they  occur  in  regions 
where  field-aligned  currents  are  high  and  where  the  auroral-electron  precipi- 
tation is  intense. 

A model  of  the  electric  potentials  in  a cross  section  normal  to  an 
auroral  arc,  which  accounts  for  most  of  the  observations  (Swift,  Ref.  33), 
is  shown  in  Figure  13.  Here,  the  magnetic  field  lines  are  taken  to  be  vertical. 
The  equipotentials  are  symmetrical  about  the  magnetic  field  line  at  the  center 
of  the  arc.  At  any  altitude,  the  potential  is  a minimum  at  the  center  and 
increases  monotonically  toward  the  sides;  hence  the  electric  field  component 
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Figure  13.  Equipotential  model  for  electric  fields  above 

auroral  ionosphere.  Lines  represent  equipotential 
surfaces;  the  magnetic  field  is  taken  to  be 
vertical. 


normal  to  the  magnetic  field  is  directed  toward  the  center  of  the  arc. 


Toward  higher  altitudes,  along  magnetic  field  lines,  the  electric  potentials 
decrease  in  magnitude;  thus,  the  electric-field  component  along  the 
magnetic  field  is  directed  upward.  Such  a field  component  tends  to  acceler- 
ate electrons  downward  and  ions  upward.  Note  that  the  electric  field  normal 
to  the  arc  increases  toward  higher  altitudes.  No  data  are  available  on  the 
closure  of  the  potential  contours  above  the  top  of  the  diagram.  It  seems 
plausible,  however,  that  the  contours  near  the  sides  continue  upward,  parallel 
to  the  magnetic  field,  and  close  in  the  conjugate  region,  maintaining  approxi- 
mate symmetry  about  the  geomagnetic  equator. 

The  work  described  below  provides  additional  information  on  the  proper- 
ties of  this  electric  field  and  the  conditions  under  which  it  develops.  The 
properties  of  the  field  are  inferred  from  recent  satellite  measurements  of 
electron  and  ion  distributions.  In  the  following  section  (4.2),  equations 
are  developed  which  describe  the  effects  of  an  electrostatic  field  along  the 
magnetic  field  on  charged  particle  fluxes.  In  section  4.3  the  observed 
fluxes  are  interpreted  in  terms  of  the  parameters  of  the  potential  model 
described  above. 

4.2  EFFECTS  OF  ELECTRIC  POTENTIALS  AND  MAGNETIC  MIRRORS  ON  CHARGED 

PARTICLE  FLUXES 

Transformations  of  particle  fluxes  will  be  described  with  reference 
to  the  diagram  shown  in  Figure  14,  which  depicts  a converging  magnetic  field 
and  an  electrostatic  potential.  Altitudes  and  magnetic  field  intensities 
are  denoted  by  z and  B,  and  electric  potentials  by  cp.  The  subscripts  s 
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denotes  the  parameters  at  the  location  of  the  satellite.  The  electric 
field  E,  which  is  directed  upward,  extends  from  eo> where  the  potential  is 
cp^fto  Zy  The  potential  at  the  satellite  is  The  subscript  t denotes 

quantities  at  the  "top"  of  the  atmosphere,  where  particle  collisions  become 
appreciable . 


The  particle  data  which  will  be  described  in  the  following  section 

were  obtained  before  and  after  the  electric  fields  were  established. 

(E) 

We  will,  therefore,  relate  the  flux,  j , when  the  field  was  present  to  the 
flux,  j,  when  the  field  was  absent. 


In  the  absence  of  E,  the  incident  flux,  at  the  top  of  the  atmosphere 
is  equal  to  the  flux  at  zr,  i.e. 


i<V  V wt>  * j<°s’  V V 

where  the  kinetic  energy  is  wfc  » 

.2  .2 

a»d  . SHL^L s . 

Bt  8s 

The  edge  of  the  loss  cone,  og  = agc>  is  obtained  from  (33)  for  cc^ 
i.e. , 

“sc  - U-'V./V*- 

If  the  backscattered  flux  at  z is  denoted  by  jD(a  , B_,  w ), 

C o C l dC 

i!/2  <.  a £ tt,  then  at  the  satellite. 


(31) 

(32) 

(33) 

tt/2  ; 


(34) 


Bs’  ws>  “ W V WBt) 

where  again  a and  or  are  related  through  (33)  in  the  ’'loss  cone," 

S L 

TJ  2 a i n - a , and  w = w_.  . 
s sc  s Bt 


(35) 
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When  an  electric  field  is  applied  along  the  magnetic  field  direction,  j 
itself  is  no  longer  an  invariant,  as  in  equation  (31).  Rather,  the  proper 
invariant  of  the  motion  is  f,  the  phase  space  distribution  function  (the 
distribution  function  used  in  Section  2).  The  two  functions  are  related  by 
the  (non-relativistic)  formula, 

j = 2 snvri.'  (36) 

It  follows  that  the  flux  at  z^  — which  equals  the  undisturbed  flux  at 


for  o.  -S  sin^B./B  ) 2 --  is  transformed  to 
i.  1 s 

w 1 

.(E)  B w')  * — j (o  , B , w ) 

J v s’  s’  s'  w J ' s c s 

(37) 

s 

wher-\ , 4-  ecp^ 

(38) 

, .2  , Ws  . 2 

and  sin  a , = — r sin  a 

t>  w S 

s 

(39) 

for  0 s a d n/2 

(40) 

(Pitch-angle  diffusion  has  been  ignored  here.  However,  it  will  eventually 
be  necessary  to  extend  the  analysis  to  cases  where  the  effects  of  diffusion 
are  of  magnitudes  comparable  to  electric  field  effects.  If  energy  diffusion 
is  negligible,  the  right-hand  side  of  equation  (3)  retains  its  form,  but  with 
the  space-like  variable  B/Bo  everywhere  replaced  by  the  transformed  variable 
. t}je  ohase  space  distribution  function  must  be  a function  of  w-eep 

B w 

o 

rather  than  w.  It  is  not  likely  that  energy  diffusion  can  be  ignored  if  the 

electric  field  is  set  up  by  plasma  turbulence  of  electrostatic-type  waves.) 

Although  the  primary  flux  for  a[  > sin  *(  is  unknown,  the  potential 

s w 

s 

cp^  will  reflect  upgoing  electrons  observed  at  with  energies  w^  and  pitch 
angles  a's  such  that, 
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, I 2 , 

W Icos  O' 
s | s 


B -B,  _ 

si  , 2 
+ — - — sin  <x 

B s 

s 


(41) 


Hence,  in  the  pitch-angle  range 

, -1  /Ws  ‘ Bs  \k  . 

sln  — • ¥7  / * * 


’ < tt/2 

s 


(42) 


(43) 


with  «p1  < wj,  ■?  av1  (•£ — jj-g-  ) 

s 1 

ar;  well  as  for  0 < Ot'  g tt/2  with  w ^ et&  , The  dov/ngoing  flux  at  7.  is  equal 

S G L S 

to  the  upfioing  flux,  i . e. 


■i(! ^ (tfa,  Bs,  w*g)  = (n  - Gd,  Bg,  w*B). 


(44) 


At  the  top  of  the  atmosphere,  zt>  through  the  potential  cpb  , the  electron 
flux  becomes 

w 


(v 

V 

wt)  * 

where 

"t  ' 

w*  + 
s 

^b 

’ ?b 

tnd 

sin2 

w* 

s 

at  * 

w1" 

+ ecp, 

t .(E) 
s 


B 


t , 2 , 

— - sin  a’  . 
B s 

s 


(45) 

(46) 

(47) 


Hence,  the  edge  of  the  loss  cone  of  the  observed  flux  at  zg  is  now  energy 
dependent  and  has  increased  from  the  value  given  in  Equation  (34)  to 
vl  + e<o.  B \ ^ 

\ 


a'  * sin 
sc 


il  K * M* 

\ w ' ' 


(48) 


The  pitch-angle  distributions  of  the  electrons  can  therefore  be  used  to 
determine  tp^. 
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Note  also  that  the  electron  flux  reflected  by  the  potential  above  the 
satellite  (Eq.  44)  is  also  transformed  to  the  top  of  the  atmosphere  accord- 
ing to  Eqs.  (45)-(47),  and  those  electrons  with 

a'  2 a'  (49) 

s sc  v ' 

are  reflected  back  toward  the  satellite  by  the  magnetic  field.  Since  Eqs.  (49) 

and  (42)  are  satisfied  by  electrons  in  the  energy  range 


ecpi  + ecpb  B1/Bc  ( ^ f__L. 

(1-B^Bj.)  " Ws  ^1  (1-B1/Bs;  » 


(50) 


such  electrons  will  be  trapped  between  the  magnetic  mirror  at  the  top  of  the 
atmosphere  (zt)  and  the  upper  boundary  (z^)  of  the  electric  potential.  The 
electron  data  discussed  in  the  next  section  indicates  that  electrons  are  indeed 
trapped  in  this  manner. 


The  AURORA  code  was  used  to  compute  the  flux  of  electrons  backseat tered 
from  the  atmosphere.  An  iteration  procedure  was  used  to  determine  the  appro- 
priate incident  flux  which  included  not  only  that  given  by  Eq.  (45)  but  also 

the  backscattered  electrons  which  were  reflected  back  to  the  atmosphere  by 

(E) 

the  electric  field.  The  latter  flux,  Jv^  , was  determined  from  the  equations 


(<*t.  Bt,  v)  » (n  - V Bfc,  w) 


where 

0 £ c*t  £ tt/2 

for  w s e92 

and 

sin 

w - eq,2 

B, 

•17) 

w 

Here 

*• 

the  backscattered  flux. 

Ba 


(51) 

(52) 

(53) 
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The  backscattered  flux  which  reaches  the  satellite  is 


w- 


Jb  4(E), 


j ; (as-  v w - ^b>  * jt«v  v «> 


(54) 


2 w 

where  sin  a'  = 

s w - e<pb 


Bs  . 2 
— sin  at 


(55) 


w > «pb  . 


4.3  ELECTRIC  FIELDS  INFERRED  FROM  SATELLITE  DATA 
4.3.1  Satellite  Measurements 

Exceptionally  valuable  data  on  particle  distributions  in  the  magneto- 
sphere are  being  obtained  with  the  satellite  1976-653,  a high-altitude  polar 
orbiting  satellite  that  has  been  launched  recently,  Drs.  E.G.  Shelley,  R.D. 
Sharp,  and  R.G . Johnson  of  LMSC  are  the  principal  investigators  of  ? major 
experiment  on-board  which  measures  the  energy  and  angular  distributions  of 
ions  in  the  energy  range  0.5  to  16  keV  and  electrons  in  the  range  70  eV-  24  ke'r. 
The  ions  are  detected  with  three  mass  spectrometers  which  sample  the  mass-per- 
unit  charge  (m/q)  distribution,  in  the  m/q  range  1-30,  once  per  second. 

Since  particle  fluxes  change  rapidly  in  space  and  time  at  high  latitudes,  it 
is  important  to  emphasize  the  time  periods  over  which  the  energy  and  pitch- 
angle  ranges  are  sampled.  The  energy-per-unit  charge  range  0.5  to  16  keV  is 
covered  by  each  spectrometer  in  12  steps  (0.5,  0.68,  0.94,  1.28,  1.76,  2.4, 

3.3,  4.5,  6.2,  8.5,  11.6  and  16.0  keV) . At  any  particular  time,  the  energy 
settings  of  the  detectors  are  at  different  values.  They  are  cycled  through  4 
values  every  64  seconds,  remaining  on  each  step  for  16  seconds.  Thus,  a 
12 -point  energy  spectrum  is  acquired  from  the  3 spectrometers  every  64  seconds. 
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The  pitch-angles  of  the  particles  are  determined  by  correlating  the 
view  directions  of  the  detectors,  which  are  perpendicular  to  the  spin  axis 
of  the  satellite,  with  the  direction  of  the  magnetic  field.  The  satellite 
is  spinning  at  3 RPM  with  its  spin  axis  perpendicular  to  the  orbital  plane. 
Thus,  essentially  the  entire  pitch-angle  range  is  sampled  by  each  detector 
every  10  seconds. 

The  electrons  were  detected  by  four  spectrometers,  designated  CME  A, 

CMC  B,  CME  C,  and  CME  D.  These  spectrometer  channels  detected  electrons 
principally  in  the  ranges  .074-. 238  keV,  .352-1.13  keV,  1.58-5.04  keV, 
and  7.3-23.5  keV,  respectively. 

*4”  -j- 

Data  from  the  first  few  orbits  revealed  several  events  of  H and  0 
ions  streaming  upward  along  the  magnetic  field  in  the  auroral  and  polar  re- 
gions, even  during  times  of  magnetic  quiet.  In  the  following  section,  we  offer 
an  interpretation  of  the  data  obtained  on  two  of  these  events  which  were  re- 
ported recently  by  Shellay  et  al.  (Ref.  5).  We  also  analyze  an  additional 
event  for  which  we  utilize  the  electron  data  as  well  as  the  ion  data  to 
interpret  the  results. 

4.3.2  Interpretation  of  Data 

In  the  two  examples  of  the  satellite  data  reported  by  Shelley  et  al. 

f + 8 2 

(Ref.  5),  approximately  equal  fluxes  of  0 and  H ions  (about  10  /cm  .sec.sr. 
keV  at  1.23  keV)  were  observed  simultaneously.  (The  detector 
settings  were  at  0.94  and  1.28  keV  at  the  time  of  the  measure- 
ments). Moreover,  the  pitch  angle  distributions  of 
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both  ions  were  nearly  identical.  They  were  sharply  peaked  along  the 
field  with  a half -width-at-half-heigbt  of  roughly  15°.  Both  events  were 
observed  at  high  magnetic  latitudes  (?*  69°);  one  was  observed  at  a local  time 
of  about  14.5  hr  at  an  altitude  of  about  7680  km,  the  other  was  observed 
at  a local  time  of  about  21hr  at  an  altitude  of  about  5450  km. 


These  data  indicate  that  the  accelerating  potentials  below  the 
satellite  were  about  1 kV.  From  the  altitudes  of  the  observations,  we  note 
that  the  electric  field  intensities  must  have  been  higher  than  about 
1000V/6000  km  or  0.2  mV/m.  In  Such  a field  ambient  H+  and  0+  would  be 
accelerated  to  1 keV  in  about  22  sec  and  90  sec,  respectively.  Now,  if 
we  assume  steady-state  conditions  on  a time  scale  of  90  sec,  the  observed 
ion  flux,  should,  as  discussed  tn  Section  4.2,  be  a factor  wi/wth  higher 
than  the  flux  of  ambient,  ionospheric  ions  at  the  base  of  the  acceleration 
region.  Here,  wfc^  sa  0.1  eV  (the  ion  thermal  energy  in  the  ionosphere),  and 
wi  ” wth  + ^ keV<  Moreover,  the  pitch  angles  of  the  energetic  ions,  Or  , 
should  be  related  to  the  pitch  angles  of  the  thermal  ions,  by  the 

equation 


0 W.,  B n 

. 2 „ th  s .2 

sin  Of,  » — sin  a., 

i w,  B th 

i o 


(56) 


where  Bg/Bo  is  the  ratio  of  the  magnetic -fie Id  intensity  at  the  satellite  to 
that  at  the  base  of  the  acceleration  region.  According  to  this  equation 
the  pitch-angle  distribution  of  the  ions  should  have  been  ever,  more  sharply 
peaked  along  the  field  than  indicated  by  the  data.  Such  a broadening  of 
the  distribution  could  have  been  due  to  scattering  of  the  ions  by  the  waves 
in  the  anomalous  resistivity  region.  The  magnitude  of  the  observed  flux. 


71 


however,  indicates  that  the  ions,  even  at  low  energies,  are  not  severely 
restrained  by  the  waves.  If  we  equate  the  ion  flow  (number/sec)  in  the 
magnetic  tube  sections  at  the  source  and  at  the  observation  altitudes, 
we  get. 


1.25 


15 


< fvi  A ~ 2r'k  dw.  j ta". , B , w.)  sin  a-,  cos  a.  da.  (57) 
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where  f is  the  velocity  distribution  function  of  the  ions,  assumed  to  be 

Maxwellian,  and  A and  A are  the  normal  cross  sectional  areas  of  the  tube 

o s 

at  the  source  and  at  B . Now,  if  we  assume  that  all  thermal  ions  with 

s 


V|  j > 0 along  the  field  become  accelerated,  then, 
*'fV|  i>  = 
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(58) 


Here,  n is  the  number  density  of  the  ions  at  the  source  and  v is 
thermal  velocity  of  the  ions.  By  substituting  this  result  and  the 
flux  in  Eq.  (57),  we  find  that 


2 x LO7 


the 

observed 


(59) 


where  we  have  put  A /A  - B /B  • Thus,  since  v . « 2 x 10  cm/sec  3000°K) 
o s o s t h 

for  0+  ions,  and  using  B^/B^  3 [,■-}  (6000  4-  R^.)^/(2500  + R^)^],  we  find 

that  an  0+  ion  density  of  about  300/cm^  at  the  base  of  the  acceleration  re- 

h 

gion  can  provide  the  observed  flux,  hi nee  the  15  ion  thermal  velocity  is 
about  4 times  that  of  the  Q+  ions,  the  required  ion  concentration  is 
about  80/crri^.  The  distributions  of  ions  in  the  polar  ionosphere  are  quite 
variable.  However,  such  concentrations  at  altitudes  of  several  thousand 
kilometers  are  not  unreasonable. 
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The  observed  pitch-angle  distribution  of  the  ions  can.  also  be  accounted 
for  by  assuming  time  and/or  space  variations  of  the  potential  model  shown  in 
Figure  13.  In  that  model  note  that  ions  are  accelerated  by  a weak  electric 
field  component  along  the  magnetic  field;  but  at  high  altitudes,  they  en- 
counter a much  stronger  electric  field  which  is  perpendicular  to  the  magnetic 
field.  There,  the  ions  will  acquire  a high  E x 3 drift  velocity.  Mozer  et  al. 
(Ref,  7),  with  instrumentation  on  the  same  satellite,  has  measured  electric 
fields  perpendicular  to  the  magnetic  field  as  high  as  1 V/m.  In  such  a field 
at  the  altitude  of  the  satellite  the  ions  would  acquire  a v component  of 
about  170  km/sec,  which  is  the  velocity  of  a 2.4  keV  0 ion.  Hence,  a per- 
pendicular electric  field  of  about  .2  V/m  for  the  0+  ions,  .7  V/m  for  H+, 
can  account  for  the  15-degree  widths  of  the  ion  distributions.  Certainly 
noth  this  mechanism  and  wave-particle  interactions  must  be  affecting  the  ion 
distributions.  The  relative  importance  of  the  roles  of  these  mechanisms, 
however,  is  uncertain  at  the  present  time. 

The  analysis  of  the  additional  event  of  this  type  that  we  investigated 
indicated  that  the  electric  potential  difference  along  the  magnetic  field 
increased  with  time  and/or  space  and  that  the  field  extended  above  as  well 
as  below  the  satellite.  The  data  were  obtained  at  105*:  UT,  September  15, 

1>76,  at  dusk,  ~ 18  hr.  LT,  near  the  geog raphi cal  coordinates  101.5°  E.  Long., 
69.5°  ft.  let.,  and  at  an  altitude  of  7^80  km.  At  the  time  of  the  measure- 
ments the  satellite  was  moving  toward  higher  latitudes. 

The  data  obtained  during  three  successive  spin  periods  of  the 

+ +• 

satellite  were  investigated.  Field-aligned  H and  0 ions,  with  pitch-angle 
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distributions  similar  to  those  described  above,  were  observed  on  the  2nd  and 
3d  spin  periods  but  not  on  the  first.  Accordingly,  it  will  be  assumed  that 
the  accelerating  field  was  not  present  at  the  time  of  the  1st  spin  period. 

(It  is  most  likely  that  the  onset  of  the  electric  field  was  due  to  the  rela- 
tive motion  of  the  satellite  and  the  region  containing  the  field.)  The 
electron  distributions  were  fairly  uniform  during  several  spin  periods  prior 
to  the  observation  of  the  ions.  The  electron  flux  measured  during  the  1st 
spin  period  will,  therefore,  be  taken  to  be  the  undisturbed  flux,  j (<>,,>  Bs,  w ), 
according  to  the  notation  used  in  Section  4.2,  and  the  electron  flux  measured 

/ E*  \ 

later  will  be  the  flux  perturbed  by  the  field,  jv  (cr',  B , w1).  The  electron 
distributions  changed  rapidly  during  the  time  of  the  measurements.  Hence, 
the  pitch-angle  data  shown  below  for  the  high  counting  rate  channels,  CME  B 
and  CME  C,  are  plotted  separately  for  each  of  tne  two  pitch-angle  scans  that 
were  made  during  a particular  spin  period.  These  individual  scans  during  a 
spin  period,  N,  are  designated  by  N for  the  scan  of  pitch  angles  from  0 to 
and  N+  for  the  succeeding  scan  from  u to  0.  The  counting  rates  of  the  CME  A 
and  CME  D channels  were  low;  to  improve  the  statistical  accuracy  of  the 
data,  those  counting  rates  were  averaged  over  similar  pitch-angle  intervals 
observed  during  scans  N and  n"*”. 

The  pitch-angle  measurements  obtained  with  the  CME  A detector  are 
shown  in  Figures  15a  and  b.  The  triangles  in  Figure  15a  designate  the 
averaged  counting  rates  obtained  during  spin  1 and  the  open  circles  designate 
the  rates  obtained  during  spin  3.  Figure  15b  shows  the  counting  rates  obtained 
on  the  individual  scans  2~  and  2+.  Note  that  the  undisturbed  flux  of  these 
low-energy  electrons,  .074-, 238  keV,  is  quite  uniform,  even  in  the  loss  cones. 
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Figure  15a.  Pitch-angle  data  obtained  on  spins  1 and  3 by  the  CME  A detector 


aNooas-jnvHAiNnoD 


the  CME  A detec 


This  isotropy  is  expected  on  the  basis  of  the  pitch-angle  diffusion  results 
discussed  in  Section  2.  A diminution  of  the  flux  of  electrons  returning 
from  the  near  atmosphere  (the  near  loss  cone,  a £ 130°)  is  indicated  in  the 
perturbed  flux  on  spins  2 and  3. 

The  pitch-angle  distributions  measured  with  the  CME  B detector  on 
spins  1,  2 and  3 are  shown  in  Figures  16a,  b and  c,  respectively.  The  count- 
ing rates  from  both  scans  of  the  pitch  angle  range  for  each  spin  period  are 
shown  in  the  Figures.  Note  that  the  electrons  heading  toward  the  near  atmosphere 
are  approximately  isotropic,  except  in  scan  3 . Also  note  that  the  near 
loss  cones  become  progressively  deeper  and  wider  in  the  successive  spin  periods, 
and  that  beyond  the  edge  of  the  loss  cone  the  counting  rates  increase  in  the 
successive  spin  periods. 

The  pitch-angle  data  obtained  with  the  CME  C detector  display  charac- 
teristics similar  to  those  of  the  CME  B detector.  For  these  higher-energy 
electrons,  however,  the  far  atmosphere  loss  cone,  a ^ 30°,  is  evident  in 
the  1 scan.  These  data  are  shown  in  Figures  17a,  b and  c. 

The  pitch-angle  data  obtained  on  spins  1 and  3 by  the  CME  D detector 
are  shown  in  Figure  18.  The  statistical  accuracy  of  the  data  is  rather 
poor.  Nevertheless,  the  near  and  far  atmospheric  loss  cones  can  be  discerned. 
Moreover,  the  counting  rates  (open  circles)  on  spin  3 exhibit  maxima  near 
the  edges  of  the  loss  cones,  at  40°  and  140°. 


77 


COUNTS/HALF  SECOND 


180  160  HO  120  100  80  60  40  20  0 

PITCH  ANGLE  (DEG) 

Figure  17a.  Pitch' angle  data  obtained  on  spin  1 by  the  CME  C detector. 
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Figure  17c,  Pitch-angle  data  obtained  on  spin  3 by  the  CME  C detector. 
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Figure  18  Pitch-angle  data  obtained  on  spins 


The  energy  spectra  of  the  electrons  detected  on  spins  1 and  3 are  shown 
in  Figure  19.  These  spectra  were  obtained  simply  by  dividing  the  mean 
values  of  the  counting  rates,  outside  the  loss  cones,  by  the  geometric  fac- 
tors of  the  detectors,  and  plotting  the  results  at  the  peak-response 
energies  of  the  detectors.  The  triangles  shown  in  the  figure  will  be  ex- 
plained later. 

The  potential  difference,  cpb>  below  the  satellite  can  be  estimated 
from  the  widths  of  the  near  loss  cones.  The  relationship  of  the  edge  of  a 
loss  cone  to  the  electron  energy  and  is  given  in  Equation  (48).  In 
Figure  20  the  width  of  the  loss  cone  given  by  this  equation  is  plotted  as 
a function  of  eqj  for  each  of  the  central  energies  of  the  electron  detectors. 
Here,  the  ordinate,  Aar'  , denotes  the  difference  180°-  o'  . In  evaluating 
Equation  (48),  the  ratio  Bg/Bt  was  set  equal  to  the  third  power  of  the  ratio 
of  the  geocentric  distance  to  the  "top"  of  the  atmosphere  (altitude  rs  300  km) 
to  the  geocentric  distance  to  the  satellite..  The  widLhs  of  the  loss  cones 
estimated  from  the  detector  data  taken  on  spins  1,  2 and  3 are  listed  in 
Table  IV.  The  values  of  ecpb  corresponding  to  these  widths,  as  determined 
from  Figure  20,  are  listed  in  Table  V. 

The  results  indicate  that  the  potential  below  the  satellite  was  0 on 
spin  1,  and  increased  to  ~ 1.5  kV  on  spin  2 and  to  ~ 3.8  kV  on  spin  3.  The 
values  of  eqsb  inferred  from  the  CME  A data,  especially  on  spin  3,  are  not 
reliable  because  of  the  extreme  sensitivity  to  the  electron  energy  of  the 
curve  shown  in  Figure  20  for  this  detector.  Moreover,  as  discussed  below. 
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FtR”re  19.  Energy  spectra  of  electrons  measured  on  spins  1 and  3.  The 
triangles  are  the  spectral  intensities  of  the  disturbed  flux 
computed  from  the  spin  1 data  assuming  a .5  kV  potential 
above  the  satellite.  A Maxwellian  distribution  with  T = .8  keV 
is  shown  for  comparison. 
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Figure  20.  Theoretical  relationship  of  loss-cone  width  observed  by  electron 
detectors  to  the  electric  potential  below  the  satellite. 
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TABLE  IV 


Observed  Values  of  Loss  Cone  Widths  in  Degrees 


Number 

CME  A 

CME  B 

CME  C 

1 

0 

21 

21 

2 

70 

39 

27 

3 

75 

60 

35 

TABLE  V 


Values  of  e<pb  in  keV  Inferred  From  Loss  Cone  Widths 


Spin  Number 

CME  A 

CME  B 

CME  C 

1 

- 

0 

0 

2 

1.2? 

1.5 

1.5 

3 

1.2? 

3.6 

3.9 

CME  D 


~ 24 


- 27 


CME  D 
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the  reduction  of  the  counting  rate  of  this  detector  on  spin  ,'l  indicates 
that  the  electrons  above  the  satellite  were  accelerated  by  ar  additional 
potential  of  about  0.5  kV.  The  detector  therefore  may  have  been  counting 
electrons  with  energies  far  from  its  principal  energy-response  band. 

The  potential  below  the  satellite  was  also  estimated  from  the  depths 
of  the  loss  cones,  by  using  the  AURORA  code  to  compute  the  flux  of  back- 
scattered  electrons  for  various  assumed  values  of  the  potential.  The  inci- 
dent flux  at  the  top  of  the  potential  region  was  taken  to  be  isotropic  and 
to  have  the  spectra  shown  in  Figure  19.  The  spin  1 spectrum  was  used  for 
the  <pb  * 0 condition,  and  the  spin  3 spectrum  was  used  for  cpb  0.  Results 
were  obtained  for  cpb  = 1.5,  3.0  and  4.5  kV.  In  general,  it  was  found  that 
the  backscattered  flux  of  electrons  which  reached  the  satellite  decreased 
rapidly  at  <pb  increased.  However,  for  the  = 0 case,  the  computed  back- 
scattered  flux  was  less  than  that  observed  at  the  satellite  on  spin  1.  This 
discrepancy  is  believed  to  exist  because  the  high-energy  AURORA  code,  which 
was  used  for  this  analysis,  does  not  compute  the  production  of  secondary 
electrons.  These  electrons  are  not  important  for  most  applications  of  the 
code,  but  they  are  significant  in  the  computation  of  the  backscattered  flux. 
Because  of  this  difficulty,  the  loss-cone  depth  computed  for  each  value  of 
<pb  was  divided  by  the  depth  for  S>b=  0 and  the  ratios  were  compared  with  the 
ratios  formed  by  dividing  the  loss-cone  depths  observed  on  spins  2 and  3 by  the 
depth  observed  on  spin  1.  This  procedure  tends  to  cancel  errors  in  the  back- 
scattered  flux  which  are  proportional  to  the  flux.  The  observed  ratios  are 
listed  in  Table  VT.  Comparison  with  the  computed  ratios  in  Table  VII  indicates 
that  tpb  « 1.5  kV  on  spin  2 and  cp^  «s  4 kV  on  spin  3.  These  values  are  con- 
sidered to  be  in  good  agreement  with  those  listed  in  Table  V. 
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TABLE  VI 


i 

■t 


Observed  Ratios  of  Loss  Cone  Depths 


Spin  Number 


CME  C Data 

.Depth  on  numbered  spin. 
Depth  on  Spin  1 


1 


1.0 


2 3.0 

3 8.8 


TABLE  VII 

Computed  Ratios  of  Loss  Cone  Depths 


AURORA  Code  Results 
Potential,  <p.,  in  kV  .Depth  for  gp  * cpK) 

’‘Depth  for  q>=0  ' 


0 1 

1.5  2.8 


H 


3.0 


5.6 


Hie  0+  and  H+  ion  data  indicate  somewhat  higher  values  of  qj^.  On  spin 

1,  of  course,  the  ions  were  not  present,  which  is  consistent  with  co^  = 0.  Hie 

data  from  the  energy-mode  channels  prevailing  at  the  time  of  spin  2 indicated 

high  fluxes  of  0+  at  3.3  keV  and  H+  ions  at  0.94  keV.  On  spin  3 high  fluxes 
+ + 

of  both  0 and  H ions  were  observed  at  4.5  keV. 

Tlie  electron  data  also  indicate  that  a potential  of  about  0.5  kV  existed 
above  the  satellite  on  spin  3.  Application  of  Equation  (37),  taking  as  the 
undisturbed  flux,  j(ao,  w^),  the  flux  detected  on  spin  1 (see  Figure  19)  and 
assuming  cp^  ss  . 5 kV,  the  calculated  flux  disturbed  by  the  electric  field  re- 
produces fairly  well  the  flux  observed  on  spin  3.  The  flux  values  calculated  in 
this  manner,  using  the  CME  A,  CME  B,  and  CME  C data  on  spin  1,  are  given  by  the 
triangles  in  Figure  19 • This  potential  is  also  consistent  with  the  peaking  of 
the  OFF  A counting  rates  at  90°  and  the  approximate  symmetry  of  the  counting 
rates  about  90°  on  scans  2+  and  3.  This  detector  is  evidently  couni  ing  principally 
electrons  which  are  backccattered  from  the  atmosphere  and  reflected  from  the 
potential  above  the  satellite.  Furthermore,  the  potential  above  the  satellite 
is  inferred  from  the  enhanced  counting  rates  of  the  CME  C and  CME  D detectors 
near  the  edges  of  the  loss  cone.  For  small  values  of  B /B-, , the  electric 
potential  above  the  satellite  reduces  the  pitch  angles  of  the  electrons  such  that 
the  lower  energy  electrons  within  the  detector  energy  bands  are  more  nearly 
aligned  with  the  magnetic  field.  This  effect  appears  tr,  account  entirely  for  the 
CME  D counting  rates.  The  peak  in  the  CME  C detector  may  also  be  due  to  trap- 
ping of  the  electrons  between  the  electric  field  above  the  satellite  and  the 
magnetic  field  below.  Application  of  Equation  (50),  with  Bj  set  equal  to  the 
magnetic  field  intensity  1000  km  above  the  satellite,  indicates  that  electrons 
of  energies  .9^  - 2.63  keV  (within  the  response  of  the  CME  C detector)  may  be 
trapped  in  this  manner  for  = 3.7  kV  and  = 0.5  kV. 
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IONOSPHERIC  IRREGULARITIES  DUE  TO  ATMOSPHERIC  HEAVE 


A high  altitude  nuclear  detonation  heats  the  atmosphere  principally 
through  the  absorption  of  X-rays  near  80  km.  The  heated  atmosphere  rises, 
approximately  vertically,  and  through  collisions  moves  plasma  across  the 
magnetic  field.  As  in  the  case  of  the  dynamo,  an  electric  field  is  induced 
in  the  plasma.  The  polarization  field  resulting  from  the  induced  field 
propagates  up  magnetic  field  lines  and  causes  the  plasma  all  along  the  field 
line  to  convect  outward  at  the  E x B drift  velocity.  The  situatioi  is  similar 
to  the  atmospheric  dynamo  which  is  driven  by  the  tidal  motions  of  the  atmo- 
sphere. If  the  velocity  of  the  ions  relative  to  the  electrons  exceeds  the 
thermal  velocity  of  the  ions,  the  Farley-Buneman  instability  (Refs,  8 and 
51)  may  occur  and  lead  to  the  formation  of  magnetic -fie Id  aligned  plasma 
irregularities. 


The  Farley-Buneman  instability  is  based  on  the  familiar  plasma  ion 
wave  instability  known  as  the  two-stream  instability.  This  instability  for 
the  case  of  a highly  ionized  collisionless  plasma  has  received  considerable 
theoretical  study  and  is  well  understood  (Refs.  51,  52,  53,  54  and  55). 

Farley  (Ref.  8)  extended  the  simple  two-stream  theory  to  the  problem  of 
field-aligned  irregularities  in  the  ionosphere,  taking  into  account  the  presence 
of  the  geomagnetic  field  and  the  effect  of  collisions  between  the  charged  and 
neutral  ionospheric  particles.  The  results  show  that  plane  electrostatic  waves, 
and  associated  irregularities  in  ionization  density  along  field  lines,  will 
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appear  spontaneously  in  ionospheric  regions  where  strong  currents  flow  per- 
pendicular to  the  geomagnetic  field  lines,  providing  4 . > a^.  Here  it 

is  the  unit  wave  propagation  vector,  is  the  relative  drift  velocity  be- 
tween the  electrons  and  ions,  and  a^,  is  the  thermal  speed  of  the  ions. 


Buneraan  (Ref.  56)  has  also  considered  the  electrojet  problem,  using 
the  Navier-Stokes  fluid  equation  rather  than  the  more  complicated  Boltzmann 
equation  applied  by  Farley.  In  the  high  collision  frequency  limit  where  the 
two  approaches  may  be  compared,  Farley  found  that  his  and  Buneman's  dispersion 
equations  are  essentially  in  agreement. 

The  conditions  required  for  the  onset  of  the  Farley-Buneman  instability 
have  been  verified  for  both  the  equatorial  and  auroral  regions  (Refs.  57,  58, 
59,  60  and  61). 


Much  attention  has  already  been  given  to  processes  which  may  be  respon- 
sible for  the  striations  photographed  during  the  high  altitude  nuclear 
tests;  and,  indeed,  the  effects  of  specific  mechanisms  have  been  incorporated 
in  various  codes.  To  our  knowledge,  the  mechanism  discussed  here  has  not  yet 
been  considered  (Ref.  62).  Yet,  it  appears  to  be  quite  straightforward,  and 
it  does  account  for  a feature  not  explained  by  other  mechanisms,  namely  the 
30-40  sec  delay  time  in  the  appearance  of  the  striations  following  a nuclear 
burst.  In  the  proposed  mechanism,  this  feature  is  expected  because  the 
atmosphere  below  the  burst  does  not  move  upward  sufficiently  rapidly  at 
timis  shorter  than  the  observed  delay  time. 


93 


In  this  section  we  discuss  the  possibility  that  the  Far  ley-Buneman 
instability  might  occur  following  a nuclear  detonation. 

The  ionosphere  and  debris  near  a high  altitude  nulcear  explosion  are 
rendered  highly  conducting.  The  net  effect  is  to  enhance  the  tendency  of 
plasma  to  move  as  a single  fluid  in  which  ions  and  electrons  are  tied  firmly 
to  the  field  lines.  The  likelihood  of  plasma  instabilities  driven  by  rela- 
tive motion  between  ions  and  electrons  is  therefore  low  in  the  ionosphere 
below  the  detonation  point.  At  large  distances,  however,  the  conductivity 
remains  only  slightly  elevated;  so  the  question  remains  whether  instabilities 
could  be  excited  there. 


A simplified  model  was  usnd  to  estimate  the  relative  drift  motion  of 

the  ions  and  electrons  in  regions  of  the  ionosphere  where  the  atmospheric 

heave  is  not  prominent.  The  model  is  illustrated  in  Figure  21.  The  ion 

velocity  across  the  magnetic  field,  driven  by  the  expansion  of  the  atmosphere, 

is  assumed  to  be  uniform  within  the  circular  region  of  radius  R.  E,  . 

ina. 

is  the  resulting  v^  x B electric  field  induced  by  the  plasma  motion.  This 
field  is  generally  directed  toward  the  west  for  the  burst  occurring  in  the 
northern  hemisphere.  Outside  the  circular  region  the  heave  motion  is  assumed 
to  be  negligible.  The  field  will  therefore  produce  polarization  charges 

at  the  boundary  of  the  circular  region,  as  indicated  in  the  figure,  and  these 
charges  will  propagate  upward  along  the  local  magnetic  field  lines.  Since 
the  atmosphere  expands  on  a long  time  scale  --  tens  to  hundreds  of  seconds  — 
charges  spread  to  distances  along  the  magnetic  field  which  are  much  greater 
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than  the  transverse  extent  of  the  electrostatic  field.  Accordingly,  the 
charges  were  assumed  distributed  in  a tubular  region,  and  a two-dimensional 
solution,  in  a plane  perpendicular  to  B,  was  obtained  for  the  ensuing 
electrostatic  potential. 

In  cylindrical  coordinates,  r,  9,  the  solution  inside  the  tubular 
region  is, 

0.  = E.  , r cos  0 . (60) 

l tnd. 

The  electrostatic  field  inside  the  tube  is  therefore  equal  and  oppositely 
directed  to  Elnd  • Hence,  the  plasma  in  the  tube  will  drift  with  the 
velocity  - Eind  x B * (v^,  which  is  the  velocity  at  which  the  ions  near  the 
base  of  the  tube  are  forced  across  the  magnetic  field  by  the  neutral  atmo- 
sphere. Hence,  the  relative  velocity  of  the  ions  and  electrons  will  be  zero. 


Outside  the  tube  the  potential  solution  is, 
k R2 

0 * Eind.  ITT”  C°S  9 

o 


(61) 


where  and  kQ  are  the  dielectric  constants  of  the  plasma  inside  and  outside 

the  tubular  region.  The  electric-field  components  outside  the  region  are, 

therefore,  - 

k4R 

(62) 


Er  - E 


kjR 

ind.  kT?  cos  8 
o 


kjR 

and  Ee*Eind.  1 T?r  8in9’ 

o 


(63) 


The  distribution  of  this  electric  field  is  depicted  in  Figure  21. 
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According  to  the  MICE  code  results  for  a Standard  SPARTAN  (200  km), 

the  ion  velocity  in  the  ionospheric  region  beneath  the  burst  where  the  ions 

are  coupled  to  the  neutrals  is  about  5 x 105  cm/ sec.  Similar  velocities  are 

given  by  the  MRHYDE  code  for  the  Check  Mate  test.  If  the  inclination  of  the 

magnetic  field  is  taken  to  be  45°  and  B u . 5 gauss,  then  the  induced  field 

is  very  high,  about  1.8  V/m.  The  dielectric  constants  are  approximately 

proportional  to  the  plasma  density.  In  the  undisturbed  E-region  of  the 

ionosphere  the  ionization  density  during  the  daytime  at  mid-latitudes  is 
5 3 

about  3 x 10  ions /cm  . Results  of  the  MRHYDE  and  MICE  codes  indicate 
E-region  ionization  enhancements  by  a factor  of  100  to  1000.  Foir  the  Standard 
SPARTAN  (200  km)  the  effective  value  of  k^/ko  is  in  the  range  10  to  100  for 
R sa  1000  km.  This  value  of  R seems  appropriate  because  the  Spartan  data 
indicated  that  there  the  ion  velocities  are  down  by  more  than  a factor  of  10 
and  the  ion  temperatures  are  nearly  equal  to  the  undisturbed  E-region  ion 
temperatures.  Us;  .ng  these  values  in  Equation  (63)  to  evaluate  E,  the  result- 
ing E x B drift  velocity,  for  B * .5  gauss,  at  r «*  2R  is  found  to  be  greater 

7 

than  about  10  cm/sec.  This  velocity,  which  is  nearly  equal  to  the  relative 
ion-electron  velocity  in  the  E-region  of  the  ionosphere, is  much  higher 
than  the  ion  thermal  velocity  («  105  cm/sec)  in  the  undisturbed  E-region  of 
the  ionosphere.  Hence,  the  Farley-Buneman  instability  criterion  is  satisfied. 

In  this  analysis  the  effect  of  the  ionospheric  conductivity  on  the 
magnitude  of  the  electrostatic  field  in  the  E-region  has  not  been  taken  into 
account.  The  conductivity  would  reduce  the  field.  However,  the  effect  is 
not  expected  to  be  severe.  In  the  auroral  zones,  electric  fields  of  about 
.06  V/m  observed  above  the  ionosphere  are  found  to  produce  relative  velocities 
in  the  E-reglon  greater  than  the  ion  thermal  velocity  and  to  form  field-aligned 
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irregularities.  Such  fields  are  smaller  than  the  one  estimated  above  by 
more  than  an  order  of  magnitude. 

It  may  be  concluded  that  the  Farley-Buneman  instability  has  not  been 
ruled  out  in  the  ionosphere  near  a high  altitude  nuclear  explosion,  but, 
because  of  elevated  conductivities  there,  this  instability  mechanism  is 
not  unlikely.  In  the  distant  E-region,  the  plausibility  of  the  Far ley-Buneman 
instability  has  been  demonstrated.  However,  the  uncertainties  in  the  calcu- 
lation make  it  imperative  that  a more  detailed  calculation  be  performed. 


98 


6.  CONCLUSIONS  AND  RECOMMENDATIONS 
It  was  found  that: 

o Large  amplitude  scintillations  can  occur  at  all  times  in  mid -latitudes, 

but  show  no  strong  correlation  with  well  known  ionospheric  and  geomagnetic 
activity  indicators. 

o The  distribution  of  trapped  particles  precipitating  into  the  ionosphere 

depends  cm  diffusion  into  the  loss  cone,  backscatter  from  the  atmosphere, 
azimuthal  drift,  and  electric  fields  aligned  along  the  magnetic  field;  all 
of  which  produce  recognizable  and  distinguishable  effects. 

o Several  characteristics  of  the  electric  fields  due  to  plasma  instabilities 

in  the  upper  ionosphere,  and  the  effects  of  these  fields  on  ions  in  the  iono- 
sphere and  electrons  in  the  magnetosphere,  have  been  inferred  from  analyses 
of  the  Satellite  1973-65B  data. 

o The  Fariey-Buneman  instability  can  be  excited  by  a high  altitude  nuclear 

explosion.  The  resulting  plasma  irregularities  in  the  ionosphere  may  extend 
to  distances  of  several  thousand  kilometers  about  the  "sub-burst"  point  in 
the  ionosphere . 

It  is  recommended  that ; 

o The  WIDEBAND  probem,  and  future  programs,  be  directed  to  make  the  best 

use  of  opportunities  to  obtain  mid-latitude  scintillation  data;  and  the  theo- 
retical efforts  be  expanded  to  further  the  understanding  of  the  underyling 
mechanisms  that  lead  to  lanosphei  instabilities  - both  in  mid  and  high 
latitudes. 
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o 


Further  Satellite  1973-65B  data  on  the  energy  and  angular  distributions 


of  keV-type  ions  and  electrons  be  analyzed  to  obtain  additional  information 
on  electric  fields  produced  in  the  auroral  and  polar  regions  and  the  condi- 
tions under  which  they  occur. 

o Satellite  1973-65B  data  on  the  energy  and  angular  distributions  of 

10-100  keV  electrons  be  analyzed  to  obtain  definitive  information  on  the 
diffusion  rates  of  trapped  electrons  at  mid-latitudes,  hence  to  calibrate 
the  theoretical  model  discussed  in  Section  2. 

o Further  studies  be  conducted  to  determine  the  electric  fields  and 

currents  in  and  above  the  ionosphere  that  result  from  high-altitude  nuclear 
explosions  and  to  test  the  results  against  existing  plasma-instability 
criteria. 
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ATTN:  FCPH 

Director 

interservice  Nuclear  Weapons  School 
ATTN:  Document  Control 

Director 

Joint  Stmt.  Target  Planning  Staff,  JCJ» 
ATTN:  Doc.  Control 

Chief 

Livermore  Division,  Field  CoMnand,  DNA 
Lawrence  Livermore  Laboratory 
ATTN : FCPRL 


DEPARTMENT  OF  1) KEENS E (Cont  i nui-il) 

Di rec  tor 

Nat iunal  Security  Agency 
ATTN : TDL 

ATTN:  Technical  library 

0.ICS/J-  i 

ATTN:  WVMCCS,  Eval.  Off.  Mr.  Toma 

ATTN;  J~3,  Knvi ronmental  Services  Div. 

OJC‘7.1-5 

ATTN:  Nuc.  Div. 

DEPARTMENT  OF  THE  ARMY 

Commando  r / D 1 r t?c  cor 
Atmospheric  Sciences  Laboratory 
US  Army  Electronics  Command 
ATTN:  R.  Olsen 

ATTN;  E.  Butterfield,  DRSEL-B1.-SY-R 
ATTN:  t\  F,,  Niles,  DRSEL- BL-SY-S 

ATTN:  Robert  Rubio 

ATTN:  H.  BalUrd 

Director 

HMD  Advanced  Tech.  Center 
Huntsville  Office 

ATTN : ATC-T,  Melvin  T.  Capps 

ATTN:  ATOO,  W.  Davies 

Program  Manager 
HMD  Program  Office 

ATTN:  DACS-BMT,  John  Shea 

Commander 

BHD  Jystem  Command 

ATTN:  BDMSC-TEN,  Noah  J.  Hurst 

Dep.  Chief  of  Staff  for  Research  Dev,  & Acq, 
DeparLheut  of  the  Army 
ATTN:  DAT.A-WSZ-C 

ATTN : DAMA-CSZ-C 

Chief  of  Engineers 
Department  of  the  Army 

ATT..:  Fernand  DePercln 

Deputy  Chief  of  Staff  for  Ops.  & Plana 
Department  of  the  Army 

ATTN:  DAMO-DDC,  COL  David  W,  Einsel 

ATTN : DAMC  *SSN 

Commander 

Harry  Diamond  Laboratories 
2 cy  ATTN:  DRXDO-NP,  F.  N*  Wimcnitr. 

Director 

US  Army  Ballistic  Research  Labs. 

ATTN:  Tech.  Lib.,  Edward  Baicy 

ATTN:  John  C.  Heater 

ATTN:  J.  Helmet l 

ATTN:  M.  .(regel 
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DUY-RTML.Yi  Oi  I HI  NAVY  (Continued) 


Commander 

I'S  Armv  Klivt  ronit  S <'or.im.inil 
ATTN:  DRS1.L-RD-P 

AT1X:  DKSU.-PL-KXV,  Hum;  A.  Ikw*ke 

ATTN:  1>KSLI.-U-!K,  Edwin  T.  Hunter 

ATTN:  DRSKL-XI.-I) 

ATTN:  Stanley  F.ronruberx 

(<>rii!i.iinlr  t 

!S  Arm  Fercign  S,  lencr  Te«  h.  Center 
ATTN:  K.  Jones 

Connunde  r 

i S Army  Matter  It*  I Dev.  6 Kc.iiliiii'ss  Cnnira.ind 
AT  IN:  DKCI.Jm:,  .1,  A.  Bender 

ATTN:  DRXCD-TI. 

Cotmr..imUr 

I S Army  M I s k i I e Command 

ATTN:  DRSM1-XS.  Chief  .Scientist 

A n N : DKSMI-ARL 

Commander 

TS  Armv  Nut  I ear  Agency 

ATTN:  MONA-WF-,  J»  her be ret 

Commander 

IS  Army  S.n»i>rn  Acetify 

ATTN:  Tfilmltiii  Library 

Chief 

l‘S  Army  Rescan  h (Mine 

A 1 TN : • UI»AKI)-I,t  Hubert  Mare 
ATTN;  ( RDAKD-DLS,  Hermann  K.  Hub] 

Commander 

Whippany  Field  Of  f Ire  (BMItSC-WS-W) 
lit* j 1 Telephone  Labors tor lew 
ATTN:  John  Turner 

Commander 

White  SandH  Missile  Range 

ATTN:  STVWS-TK-NT,  Marvin  1'.  Squires 
ATTN:  AMSI-.b-NI.-sn,  Willis  W.  Webb 

DEPARTMENT  UK  THE  NAVY 

Chiel  ui  Naval  Operations 
Navy  Department 

ATTN:  OP  ytfl 

ATTN:  Code  604CJ,  Robert  Places i 

Chief  of  Naval  Re  sea  re  It 
Navy  Department 

ATTN:  Code  461 , Jacob  L.  Warner 

ATTN:  Code  461,  R.  Cracen  Joiner 

Commander 

Naval  Klee t runic  Systems  Command 
Naval  Electronic  Systems  Command  Hqs. 

ATTN:  El, EX  03 

ATTN:  PME  117-T,  Satellite  Coma.  Project  Off. 

Command  1 hr  Officer 
Naval  Intelligence  Support  Center 
ATTN:  Code  40A,  E.  Blase 

ATTN:  Document  Cant  r**l 


Commander 

Naval  Ureau  Systems  (.enter 


ATTN : 

(.'lull' 

2200. 

Jurgen  Richter 

ATTN : 

William  F. 

Mu  1 e r 

ATTN : 

Code 

2200. 

Richard  Pnppcrt 

ATTN: 

Quit' 

2200. 

I i nil  Rut  limn  1 1 cr 

TN: 

Tech. 

Lib. 

for  T.  J.  Ke.iry 

Superintendent  (Code  I -T  J /*  ) 

Naval  I'usf  graduate  School 

ATI N ; C.ule  61,  DY,  .1.  N.  Dyer 

ATTN:  Code  61,  MN , R.  A.  Milne 

ATTN:  Code  J 1 .’4 , Ter Jt.  Rpts.  Librarian 

Ml  re*  t «»r 

Naval  Research  Liberal orv 

ATTN:  Code  4004,  Kraanua  i !..  Branca  to 

A11N:  Code  5174,  Frederick  W,  Raymond 

ATTN:  Code  7l23.7,  Donald  J.  Horan 

ATT N ; Code  6623,  Francis  J.  Campbell 
ATTN:  Code  660  Jl‘,  Richard  L.  Stilt  l cr 

ATTN:  Code  2600,  Tech.  Lib. 

ATTN:  (udo  712/,  Charles  Y.  Johnson 

ATTN:  Code  71 2H,  James  I).  Kurfess 

ATTN;  W.  Neil  Johnson 
ATTN:  Code  7120,  ('.era Id  H.  Share 

ATTN:  Code  7701  , Jack  D.  ftrown 

ATTN:  Code  7709,  Wahah  All 

ATTN:  Code  7730,  Klaus  Main 

ATTN:  Code  7 730,  J.  I odder 

ATTN:  Code  7730,  J.  Davis 

ATTN:  Code  7730,  S.  1 . Ossukow 

ATTN:  Cotie  7120,  Robert  l..  K Inzer 

ATTN:  Code  7700,  Timothy  P.  Coffey 

ATTN : Code  7750,  Darrell  Strobel 

Commander 

Naval  Surface  Weapons  Center 


ATTN : 

Code  WA501 

Navy 

Nuc.  Prgms.  Off 

ATTN: 

Code  WX21 , 

Tech. 

Lib. 

ATTN: 

D.  L.  Land 

ATTN: 

L.  Rod  1 in 

Commander 

Naval  Weather  Service  Command 

Naval  Weather  Service  Headquarters 
ATTN:  Mr.  Martin 

Commanding  Officer 

Navy  Astronautics  Croup 
6 ty  ATTN:  W.  J,  Cleason 

Commanding  Officer 

Nuclear  Weapons  Training  Center  Pacific 

Naval  Air  Station,  North  Island 

ATTN:  Nuclear  Warfare  Department 

Commander 

Nuclear  Weapons  Training  Croup  Atlantic 
ATTN:  Nuclear  Warfare  Department 

DEPARTMENT  OF  THE  AIR  FORCE 

Commander 

ADC/DO 

ATTN:  DOF 

ATTN:  DOFS 

ATTN:  D0A 
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PI  FAKEMENT  OF  THE  AIK  FORCE  Ul..nt  intird) 

Commander 

ApcnM/XJ’D 

ATTN:  Xr^I>i> 

ATTN:  Hqs.  1'ith  Arrcsp.ii  »•  I*»rre  (ITS) 

ATTN:  XI\M>0.  MA ) 0.  Kurh 

AF  •■i,«»physi»-f»  .it cry , ATS l’ 


AT’I  N 

nPR, 

A) 

o.i  1.  Stair 

AIIS 

1 K|l, 

Wi 

1 i l .in  Swider , Jr  . 

a i 'i  n 

1 KB. 

A. 

1 .i  i re 

ATTN 

1 . K 

•it 

A I i N 

FUG, 

.1 . 

C.  Me t; lay 

A 1 TS 

I. KB. 

Kc 

nneth  S.  W.  Chant 

A TIN 

01*. 

John  S.  Caring 

ATTN 

Lid*. 

Horen  S.  Narclsi 

ATI :% 

LKK  , 

F. 

K.  limes 

ATTN 

I.KO, 

K. 

bert  K.  Hut  I man 

ATTN 

cl’R, 

y . 

DeKi'recu 

ATT  N 

OPR, 

J. 

Kennenly 

A T TS 

opR , 

T. 

Connol ly 

ATTN 

sroi.A , 

1..  L.  Cunh.i 

ATTN 

i U‘K- 

1 , 

J.  C.  Cl wick 

ATTN 

i.KB, 

K. 

K.  Innes 

ATTN 

LK«, 

T. 

Kerne she a 

ATTN 

OPR, 

F. 

Hill  IngL'ly 

AITS 

l.KS, 

K. 

A.  Vanias se 1 

ATTN 

OPR, 

k. 

E.  Murphy 

ATT  N 

OPR. 

H. 

Gardener 

ATTN 

(»PR. 

R. 

O'Neil! 

AE  Weapons  Laboratory,  AFSC 
AITS:  SIT. 

ATTN:  I)VTt  Capt  Cary  D.  Cable 
ATTN:  DVT,  (.apt  lie  liars 
ATTN:  CA,  Arthur  H.  Cuvnt her 

ATTN:  DYC,  Joseph  Jatmi 

ATTN:  t:C,  Co!  G.  J.  Frey  or 

Al-TAC. 

ATTN:  Teeli,  hlb. 

ATTN:  TF/MaJ  Wiley 

ATTN : TN 

ATTN:  TN-E,  I.t  Col  E.  C.  Thomas 

Head qua r tern 

Air  Force  Systems  Command 
ATTN : DCS 

ATTN:  1.1  Col  B.  C.  Nealy 

ATTN;  Lt  Col  W.  Bear 
ATTN:  Ct  Col  H.  Simmons 

ATTN:  UUCP,  Capt  Koine! 

ATTN:  SCTSif,  Wpns.  b Wpns.  Effects  Div. 

ATTN:  SDS/DL5E 

ATTN:  DC1.AW 

head quarters 

Electronic  Systems  Division,  (AFSC) 

ATTN:  James  Whelan 

ATTN:  W,  C.  Horton 

Headquarters 

Electronic  Systems  Division/ ES 
ATTN:  fcSTT 

Commander 

Foteign  Technology  Division,  AFSC 
ATTN:  HI  CD.  Library 

ATTN : TDPSS*  Kenneth  N . Williams 


DEPARTMENT  OF  THF  AIR  FORCE  (Continued) 


Hq.  ESAF/IN 

ATTN : 

IN 

llq.  TSAF/RD 

ATTN : 

RDPH 

ATTN : 

RDSA 

ATTN : 

RDPS 

ATTN : 

RD‘jSM 

ATTN : 

RD 

Commander 

Hume  Air  Devi 

lopment  Cent 

it,  AFSC 

ATTN : 

OCSK,  J.  J. 

Simons 

ATTN: 

KMTI.D,  Dor. 

Library 

ATTN: 

V.  Covne 

SAMSO/AW 

lh  iy  ATTN:  AWW 

SAMSO/SZ 

ATTN:  SZJ , Major  Lawrence  Doan 


Conmiander  In 

Chief 

Strategic  Ai 

r Command 

ATTN : 

XPFS,  Ma)  Brian  G. 

Stephan 

ATTN: 

ADW,  Lt  Col  .1 . S. 

Greene,  Jr 

ATTN : 

ADWAT E , Ca pt  iirueo 

Bauer 

rSAM-.TAC/CB 

ATTN : 

CBI.  Mr.  Creasi 

Commander 

ASD 

ATTN : ASD-YH-EX,  Lt  Col  Robert  Leveret te 

DET  U 12  KS 

Space  Forecasting  Section 

ATTN:  Capt  U.  M.  lies]  in 

EN_KK( ; V ^STjAKCh  C.  DEVELOPMENT  ADMINISTRATION 

Division  of  Military  Application 
CS  Energy  Research  t,  Dev.  Admin. 


ATTN : 

Doe . 

Con . 

for 

Frank  A. 

Ross 

ATTN : 

Doe . 

Con. 

for 

Co  1 one  1 T 

. cross 

ATTN : 

Doc . 

Con. 

for 

David  H. 

Slade 

ATTN : 

Doc  . 

Con. 

I or 

Major  D. 

A.  Haycock 

I'nivraity  of  California 
Lawrence  Livermore  Laboratory 

ATTN:  Tech.  Info.  Dept.  L“ 3 

ATTN:  .!.  F.  Tinney,  L-517 

Sandia  Laboratories 
Livermore  Laboratory 

ATTN:  Doe . Con.  for  Thomas  H.  Cook,  Org.  8000 

Sand i h Laboratories 


ATTN : 

Doe . 

Con. 

for 

W.  D.  Brcwn,  Org.  1353 

ATTN : 

Doe . 

Con . 

for 

J.  C,  Eckhardt , Org*  1250 

ATTN: 

Doe . 
52  30 

Con . 

for 

Clarence  R.  Mehl,  Org. 

ATTN: 

Doc  • 

Con . 

for 

Org.  9220 

ATTN: 

Doc. 

5720 

Con. 

for 

Morgan  L.  Krasa,  Org. 

ATTN: 

Doe. 

Con . 

for 

3422,  Sandia  Rpt . Coll. 

ENERGY  RESEARCH  A DEVELOPMENT  ADMINISTRATION 
( tojn  tinui»d> 

Los  Alamos  Sclent  i f ic  Laboratory 


ATTN: 

Doc . 

Gnu  • 

for  J.  Judd 

ATTN : 

Dnc , 

Con . 

for  K.  A.  Jeffries 

ATTN: 

Doc , 

Ton. 

for  Martin  Tierney,  J-10 

ATTN : 

Dor  . 

Con . 

for  E.  A.  Bryant,  <NO  DAS I AC) 

ATTN : 

Doc  . 

Con. 

lor  Donald  Kerr 

ATTN: 

Dor  - 

Con . 

for  David  W.  Steinhaus 

Al  IN: 

Doc  . 

('nil . 

for  K.  W.  Hones,  Jr. 

ATTN : 

Rot  or cn cl* 

Library,  Ann  Beyer 

ATTN : 

pnr . 

Con . 

for  S.  Rockwood 

ATTN : 

Doc  . 

Coil. 

for  Wi 1 1 iam  Maier 

ATTN : 

Doc  • 

Con. 

fo*  John  7. inn 

ATTN : 

Doc . 

Con. 

for  W.  D.  Barfield 

ATTN: 

Doc  . 

' n . 

for  Marge  Johnson 

ATTN : 

Dm  . 

L>n . 

for  W.  M.  Hughes,  MS  559 

ATTN : 

Doc  . 

Con . 

for  John  S.  Malik 

ATTN : 

D.  M 

. Rohrer , MS  518 

l-S  Energy  Research  h Dev.  Attain* 
Divisl on  of  Headquarters  Services 
Library  branch  C-043 


ATTN: 

Doc . 

Con. 

tor 

Richard  J.  Kandel 

ATTN : 

Doc. 

Con . 

for 

H.  H.  Kurzweg 

ATTN: 

Doc . 

Con. 

for 

D.  Kohl  stud 

ATTN: 

Doc  . 

Con . 

for 

George  L.  Rogosa 

ATTN : 

Doc  - 

Con. 

for 

Class  Tech.  Lib. 

ATTN: 

Doc. 

Con . 

for 

Rpts. 

ATTN: 

Doc . 

Con. 

for 

Joseph  ».  La Fleur 

OTHER  ^VE-Jt>W_ENT  AGKNC  I ES 

Central  Intelligence  Agency 

ATTN;  RD/SI,  Rm.  5C48  Hq.  IJldg.  for  NED/OSI  - 
2G48  Hqs . 

Department  of  Commerce 
Notional  Bureau  of  Standards 


ATTN: 

Sec . 

Officer 

for 

M.  Krauss 

ATTN : 

Sec . 

Of  f icer 

for 

J.  Cooper 

ATTN: 

Sac. 

Officer 

for 

l).  R.  Lida 

ATTN: 

Sec . 

Of  f icer 

for 

Levin  H.  Gevantoum 

ATTN: 

Sec . 

Officer 

for 

Office  of  Director 

ATTN : 

Sec . 

Officer 

for 

M.  Sc hoc r 

ATTN : 

Sec . 

Of  f icer 

for 

D.  Garvin 

ATTN: 

Sec . 

Officer 

for 

B.  Stiener 

ATTN: 

Sec 

Of  f icer 

for 

George  A.  Sinnatt 

ATTN : 

Sec 

Of  f icer 

tor 

K.  Keesler 

Department  of  Commerce 
Office  of  Telecommunications 
Institute  for  Telecom  Science 
ATTN:  William  F.  Utlaut 


OTH ER  GOV FKNM ENT  AG I NCI ES  (Coni i nued ) 
NASA 

George  C.  Marshall  Space  Flight  Center 


ATTN : 

W.  T.  Roberts 

ATTN: 

R . Chappe] 1 

ATTN: 

K.  D.  Hudson 

ATTN: 

N.  H . Stone 

ATTN : 

C.  R.  Bale  her 

ATTN : 

W.  A.  Oran 

NASA 

ATTN: 

K.  A.  Sc-hiffer 

ATTN: 

M.  Tepper 

ATTN : 

D.  P.  Cau f t' man 

ATTN : 

G.  Sharp 

ATTN : 

J . Findlay 

ATTN: 

A.  Srhardt 

ATTN: 

Code  ST,  W.  L.  Taylor 

ATTN: 

P . r.a  t on 

ATTN : 

Rt  Pel  lows 

ATTN: 

Code  ST,  D . R.  Will iam 

ATTN: 

E.  R.  Sc  timer  ling 

ATTN: 

J.  Naugle 

NASA 

Johnson  Spat 

e Center 

ATTN: 

Owen  Gar  riot 

ATTN: 

Class  if  i ed  Lib.  * Code  BMf> 

NASA 

Ames 

Research  Center 

ATTN; 

C.  P . Sonet t 

ATTN: 

N-245-3,  Palmer  Dyal 

ATTN: 

M-254-4,  H.  Whitten 

ATTN: 

N- 2 4 5-4,  I 1 la  C,  Poppofl 

NASA  Scientific  A Technical  Info.  Fa-:. 
ATTN:  SAR/DLA- T85 

ATTN:  Acq.  Branch 

National  Oceanic  6.  Atmospheric  Admin. 
National  Weather  Service 
Department  of  Commerce 
ATTN:  L.  Machta 

ATTN:  J.  K.  Angell 

National  Oceanic  & Atmospheric  Admin. 
National  Ocean  Survey 
Depa r t men t of  Common: e 
Washington  Science  Center 

ATTN:  John  W.  Townsend,  Jr. 

ATTN:  Ronald  L.  Lavoie 

ATTN;  Edward  S.  Epstein 
ATTN:  Wayne  McGovern 

ATTN:  Cere. Id  A.  Peterson 


NASA 

Goddard  Space  Flight  Center 
ATTN:  R.  F.  Benson 

ATTN:  H.  A.  Taylor 

ATTN:  M.  Sugiura 

ATTN:  S.  J.  Bauer 

ATTN:  James  I.  Vette 

ATTN:  J.  Siry 

ATTN:  Technical  Library 

ATTN:  Gilbert  D.  Mead,  Code  641 

ATTN:  C.  Levin 

ATTN:  A.  C.  Aiken 

ATTN:  A-  Tempkin 

ATTN:  J.  P.  Hepnar 


National  Oceanic  A Atmoapneric  Admin. 
Environmental  Research  Laboratories 
Department  of  Commerce 

ATTN:  A.  H.  Snailey 

ATTN:  RAJ,  Donald  J . Williams 

ATTN;  RX1.  Robert  W.  Knecht 

ATTN:  Robert  B.  Doeker 

ATTN:  R43 , Herbert  H.  Sauer 

ATTN:  Eldon  Ferguson 

ATTN:  Wall her  N.  Spjeldvik 

ATTN:  George  C.  Reid,  Aeronomy  Lab. 
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OTtt HR  GOVERNMENT  AGENCIES  (Continued)  DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

National  Science*  Fnundat  Ion  Univ.  of  California  at  l.os  Angeles 


ATTN: 

W. 

D.  Adams 

Off  Ic 

e of  Contract  6 Grant  Admin. 

ATTN: 

ir 

ed  D.  White 

ATTN:  T.  A.  Far U-y,  Space  Set.  Ltr 

ATTN: 

M. 

K.  Wilson 

ATTN: 

Rolf  Sinclair 

Univ. 

of  California  at  Riverside 

ATTN : 

R. 

S , Zapolsky 

ATTN:  fames  N.  Pitts,  Jr. 

ATTN : 

K. 

P.  Todd 

ATTN:  Alan  C.  Lloyd 

AT  TN: 

W. 

H.  Cramer 

AT  TN : 

K. 

Manka 

Univ. 

of  California  at  San  Diego 

ATTN : 

Y. 

Gilman  Blake 

ATTN:  C.  E.  M<- 11  wain,  Physics  Dept 

ATTN:  S.  C.  Lin 


l!S  Arms  Control  & I > i surmair-ent  Agcy, 

ATTN:  Office  of  Director  University  of  California 

ATTN:  Reference  Info.  Center  Berkeley  Campus  Room  T 1 H 

Sproul  Hall 

pmRTOENTjlF  DEFENSE  CONTRACTORS  ATTN : Kinsey  A . Anderson 

Aerodyne  Research,  Inc.  State  of  California  Air  Resources  Hoard 

ATTN:  M.  Cjmac  ATTN:  Leo  Znfonte 


Agronomy  Corporation 


ATTN: 

S. 

A.  Bnwhill 

•space  Cor  pc 

• rat  ion 

ATTN: 

Julian  Reinheimer 

ATTN: 

J. 

Sorrels 

ATTN: 

G. 

Millburn 

ATTN: 

V. 

Josephson 

ATTN: 

M. 

Be'S 

ATTN: 

L, 

W.  Alike rman 

ATTN : 

V. 

. Wall 

ATTN: 

J. 

Woodford 

ATTN: 

s. 

P.  Bower 

ATTN: 

V. 

Keller 

ATTN: 

R. 

D.  Kuwcliife 

ATTN: 

George  Paulinas 

ATTN: 

T. 

Friedman 

ATTN: 

D. 

C.  Pec ka 

AT'IN: 

J. 

Stevens 

ATTN: 

J. 

R.  Blake 

ATTN: 

N. 

Cohen 

ATTN: 

Harris  Mayer 

ATTN: 

Thomas  D.  Taylor 

Analytical  Systems  Engineering  Corp. 
ATTN:  J,  A.  Caruso 


University  of  California 

Space  Sciences  Laboratory 
ATTN:  Bodo  Par.nJy 

Cal  span  Corporation 

ATTN:  Robert  A.  Fluegge 

ATTN:  M.  C.  Dunn 

ATTN:  Technical  Library 

ATTN:  W.  Wurster 

ATTN:  C.  E.  Treanor 

University  of  Colorado 

Office  of  Contracts  f.  Grants 

ATTN:  Jeffrey  B.  Pearce,  LASP 

ATTN:  Charles  A.  Barth,  LASP 

ATTN:  A.  Phelps,  JIU 

ATTN:  C,  Lineberger,  JILA 

The  Trustees  of  Columbia  University 

City  of  New  York 

ATTN:  Sec.  Officer,  H.  M.  Foley 

ATTN:  Richard  N.  Zare 

Comsat  Laboratories 

ATTN:  Document  Control 


Battellc  Memorial  Institute 

ATTN:  Donald  J.  Ha  mourn 

ATTN:  Radiation  Eff.  Info.  Center 


The  Boeing  Company 

ATTN:  Clen  Keister 


The  Trustees  of  Boston  College 
Chestnut  Hill  Campus 

ATTN:  R.  Hegblom 

ATTN:  R.  L.  Cii  i ovil  laa 

2 cy  ATTN:  chairmen,  Dept,  of  Chem. 

ATTN:  R.  Eather 

ATTN:  Library  for  Technical  Documents 


Brown  Engineering  Company,  Inc. 
ATTN:  J.  E.  Cato 

ATTN:  J.  Dobkina 

ATTN?  Roland  Patrick 
ATTN:  Romeo  A.  DeLiberis 


Epsilon  Laborator ieo,  Inc. 

ATTN:  Henry  Miranda 

ATTN:  Carl  Ac car do 
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